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INTRODUCTORY REMARK 



We found it convenient to make some minor changes in the positioning of 
the parts of all orders of Table II, and in the specific effect of two among them. 

First: We determine that in each order of Table II the memory position 
number x occupies not the 12 left-hand digits (cf. the second remark of 7.2, cf. also 
8.2), but the 12 right-hand digits. I.e. the digits 9 to 20 or 29 to 40 (from the 
left) for a left-hand or a right-hand order, respectively. 

Second: We change the orders 18, 19 (the partial substitution orders xSp, 
xSp' ) , so that they substitute these new positions, and these crosswise. I.e. 18 
replaces the 12 right-hand digits of the left-hand order located at x (i.e. the digits 
9 to 20 [from the left]) by the 12 digits 29 to 40 (from the left) in the accumulator. 
Similarly 19 replaces the 12 right-hand digits of the right-hand order located at x 
(i.e. the digits 29 to 40 [from the left]) by the 12 digits 9 to 20 (from the left) 
in the accumulator- 

For the standard form of a position mark, x^ = 2~^®x + 2'^®x, as intro- 
duced in 8.2, these changes compensate each other and have no effect. Therefore all 
the uses that we made so far of these orders are unaffected. 

On the other hand the new arrangement permits certain arithmetical uses 
of these orders, i.e. uses when the position x is not occupied by orders at all, but 
when it is used as (transient) storage for numbers. In this case the new arrangement 
provides very convenient 20-digit shifts, as well as certain other manipulations. 
These things will appear in more detail in the discussion immediately preceding the 
static coding in 10.7. 



10.0 CODING OF SOME SIMPLE ANALYTICAL PROBLEMS 



10^1 We will code in this chapter two problems which are typical 
constituents of analytical problems, and in which the approximation character o£ 
numerical procedures is particularly emphasized. This is the process of numerical 
integration and the process of interpolation , both for a tabulated one-variable 
fiinction. Both problems allow a considerable number of variants, and we will, of 
course, make no attempt to exhaust these. We will nevertheless vary the formulations 
somewhat, and discuss an actual total of six specific problems. 

10.2 We consider first the integration problem. 

We assume that an integral J f(z)dz is wanted, where the function f(z) 
is defined in the z-interval 0, 1, and assumes values within the range of our machine, 
i.e., size < 1. Actually f(z) will be supposed to be given at N+1 equidistant places 
in the interval 0, 1, i.e., the values f(h/N), h = 0, 1,...,N, are tabulated. We wish 
to evaluate the integral by numerical integration formulae with error terms of the 
order of those in Simpson's formula, i.e., 0(1/N^). 

This is a method to derive such formulae. Consider the expression 



(1) <p(u) = ^ / f(z)d2 . [2f(b)£ + 1 (f(liii) . 2f(b) + f(M) g« ] u . 

/ N 3 N N N N 

hz£u / 
N "^ 

There will be h = 0, 1,. . . ,N, g g g 1, | u I 1. 

Simple calculations, which we need not give here explicitly, give 
(p(0) = <p'(0) = (p"(0) = (indeed (p^"(0) = and for £=1 even (p""(0) = 0), implying 

(2) cp(u) = i^du, i^^du^ J^-du3cp-'(u3), 
and further 



(3) <?"'(u) = [f"4lii^) + f"(li^£Ji)] Si - [f"(iiiii) + f"(ii^)] ^ 

N N Na N N N3 



I [f".(lrbi) _ f..,(j^.,jj] ^ 
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(3) transforms into 



(4) 



>;3 

cp'"(u) = [f"(tiai) . £<<{hhi) + f'"(hiu) (i-g) u ] £_ + 

N N N N N^ 

+ [f"(hr^) . f-i(h:ll) - f'"(h:Ji) (l.g) U ] ^ 



N 



N 



N 



N N3 



N N 3 ^ N* 



Putting 



Mf"" =Max |f""(z)| , 



we see that (4) yields 



(5) lq)'"(u)| I 2Mf"" 0-. 2£) ^ , 



and by (2) 
(6) 



Uu,U<Me(i.2a) Is; 



Putting u = 1, and recalling (1), we find: 



(7) 



^ f(z)dz = [2f(^)g.i (f(h±l) - 2f(^) .f(t^)) £^]1.,, 

N t/ 

l,|.Jii^,7.6aj»i,. 



Putting £ - 1 and summing over h - if+1, Jf+3,..., k-1 (Jf, k =0, 1^..., N, 
k-/l^>0 and even) we get Simpson's formula: 

' f 



(8) 



Mf 



f (^) + 4f (^) + 2£i'^) + 4f(£i^) + . . . + 

N N N N 

+ 4f (k^) + 2f(k^) + 4f(li^) + f(k) 

_ N N N N 

i=i < M£!:: _i 



3N 1 



90 av^ 



180 N" 



Putting E, =^ and sumraing over h = I'+l, 1+2,..., k-l il, k = 0, 1,...,N, 
k-l'>0 and of any parity) we get the related formula: 






P 



,, ,, f(2)ds = 2 [f(!l*l) + 22f(S^) + i{h=l)] __1 + cp = 



N 



«/ 



^ ^ ^ 180 N^ 180 N^ 



10.3 We evaluate first by Simpson's forfflula^ i.e., by (8), and in 
order to simplify matters we put ^ ~ 0, k = N. Hence N must be even in this case 

PROBLEM 10 : 

The function f(z), z in 0, 1, is given to the extent that the values f(h/N), 
h "= 0, Ij.o.j N, are stored at N "** 1 consecutive memory locations p, p+1, . . . ^ p^N. It is 
desired to evaluate the integral J f(z)dz by means of the formula (8). -.»-- 



We could use either form of (8)'s right side for this evaluation. Using 
the first form leads to an induction over the odd integers h = 1, 3,...,N-1, but requires 
at each step the preliminary calculation of the expression 

[(f(h+l)/N) + 4f(h/N) H- f ((h-l)/N)]/3N. Using the second form leads to an induction 
over the integers h == 0, 1,..., N, and requires at each step the simpler expression 
(2/3N)[f(h/t'^) ] multiplied by a factor ^ = /4 or 1 or 2. This factor is best handled by 
variable remote connections. A detailed comparison shows that the first method is about 
20% more efficient, both regarding the space required by the coded sequence and the time 
consumed by the actual calculation. We will nevertheless code this Problem according to 
the second method, because this offers a good opportunity to exemplify the use of vari- 
able remote connections. 



Thus the expression to be computed is 

h=o ^ m N 



= 'A for h = 0, N 
e. I = ll for h ;^ 0, N and f h even 
[=2j -[hodd 
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This amounts to the inductive definition 



J,i = o» 



•^h =Jh.l^\ |j f(^) forh=0, 1,...,N, 
J = % 

Let A be the storage area corresponding to the interval of locations 
from p to p + N. In this way its positions will be A. 0, 1, . . , ,N, where A.h corres- 
ponds to p ■*'ii, and stores f(h/N), These positions, however, are supposed to be 
parts of some other routine, already in the machine, and therefore they need not 
appear when we pass (at the end of the whole procedure) to the final enumeration 
of the coded sequence. (Cf. the analogous situation in Problem 3.) 

Let B be the storage area which holds the given data (the constants) 
of the problem p, N. (It will actually be convenient to store p, N as p^, 
(p-^N-l) )• Storage will also have to be provided for various other quantities 
(0, 1 ; it is convenient to store 2/^; the exit- locations of the variable remote 
connections), these too will be accomodated in B. Next, the induction index h 
must be stored, h is really relevant in the combination p + h, which is a position 
mark (for A.h, storing f(h/N)), so we will store (pr^'h)^ instead. (Cf. the analo- 
gous situation in Problem 3.) This will be stored at C. Finally the quantities 
which are processed during each inductive step will be stored in the storage area 
D. 

These things being understood, we can now draw the flov diagram, as 
shown in Figure 10.1. It will prove to be convenient to store (2/3N) f(h/lM) 
after II as well as \ (2/3N) ilh/N) after III and IV (but not after V) not at D.2, 
but in the accumulator, and to transfer it to D.2 in VI only (but also in V). llie 
disposal of a^, Og, as, shovn in IX, VIII, VI will be delayed until X, and they 
will be held over m the accumulator. Finally, the transfer of (A.h+1)^ into C is 
better delayed from X until after the Cc order in XI, whereby it takes place only 
on the - branch issuing from XI, but this is the only branch on which it is needed. 
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(I> 



1 



O toD.l 



(sH° 



# 



k= N 



5i.+ -A^ 



(fc-N). 



C CA.k-fil 

— p__l 



31 



G M.0), 
D.1 



2 



0-*-h 



1 



3,3 



3n 



h-*-K*l 



aC" (Kj i-f k = 0,N 

<3L * OL^ if k ^ 0,N and ei/ert 

OC = <Ky if h odd 



^'4 



377 



%■ f(-^) *o D.Z 



TO. 



51 



o( =0C, 



fi.r 






6y 



6^f (^; to g.^ 



^f ('■^; to j>.i. 



SS. "'r 



•*^*—D.Z£h4,f(§) 



Ck-CN-l)J 



6 



-v7 7f + 



3 



8 



jx:| 



— ^'^^f(^) 



© 



5 

-*♦- 






-© 



■o 



-< («: 



A.h 



^r^) 



A 



Figure 10.1 
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B.l 

I.l 

2 

B.2 

1.3 

4 

B.3 

1.5 

6 



II. I 

2 
3 



B.4 
11,4 

5 



(to II, I) 



The static coding of the boxes I-XI follows: 



Ac 
D.l 

Ac 
C 

Ac 
11,5 

Ac 
11,3 

R 



[ 



III.l 








B.1 




D.l 


s 


Po 




B.2 




C 


s 


(ai)o 




B.3 




11,5 


Sp 


C 




II. 3 


Sp 


- 


R 


p-4i 


R ] 


2 




h 




B.4 


X 


. 


C 


a 


c ] 




R 



IV 



(to VI, 1) 
(to VI,1) 



v,i 




L 


2 


D.2 


S 


(to VII, 1) 
VI.l 


D.2 


S 


B.5 
Vl,2 


(03)0 
B.5 




(toX,l) 
VII, 1 


C 




B.6 
VII, 2 


(p^N-l)o 
B.6 


h- 


3 
(to VIII, 1) 
B.7 


IX, 1 
(02)0 


Cc 


VIII, 1 


B.7 





IX, 1 



(toX.l) 
L 
(toX.l) 



B.3 



Ac 
Ac 

Ac 
D.2 

D.2 

Ac 

Ac 

Ac 

Ac 
Ac 






Po 
Po 

(ai)o 

(p+h)^ 
P + h 

f(ii) 

N 



3N N 



^ 3N N 3N N 



£ 1- f(h)= 4_f(h) 

'^ 3N N 3N N 

e 2_ f(h) 

'^ a^ N 

£ 2_ f(h) 

^ 3N N 



(013)0 
<p+h), 

(h-(N-l))^ 

((12)0 
(ai)o 
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X,l 

2 


11,5 
D.l 


Sp 


3 


D.2 


h 


4 


D.l 


S 


(toXI.l) 
XI, 1 


C 




2 


B.6 


h- 


B.8 
XI, 3 


B.8 


h- 


4 
5 


e 
C 


Cc 


6 


B.8 


h 


7 


C 


S 



(to 11,1) 



II.5 
Ac 

Ac 
D.l 

Ac 
Ac 

Ac 

Ac 
Ac 
C 



a c 

J = I +6 -^ f(— ) 

h h-1 h 3f^ [^ 



(p+h)„ 
(h-(N-l))^ 

(h-N)„ 

(p4h)^ 
(p+h +1)^ 
(p4h-H)° 



Note, that box IV required no coding, hence its irranediate successor (¥1,1) 
must follow directly upon its immediate predecessor. However, this box has actually no 
immediate predecessor; IV, 1 corresponds to otg- Hence it must be replaced there by VI, 1. 

The ordering of the boxes is I, II; III, VI, X, XI; V, VII, VIII; IX (IV 
omitted, cf . above) and X, X, II must also be the immediate successors of VIII, IX, XI, 
respectively. Ihis necessitates the extra orders 



:ii.2 


X.l 


c 


IX,2 


X,l 


c 


XI, 8 


11,1 


c 



(Xi, (X2, OLq, correspond to 111,1, VI, 1 (instead of IV, 1), V, 1. In the 
final enumeration the three a's must obtain numbers of the same parity. This may 
necessitate the insertion of dummy (ineffective, irrelevant) orders in appropriate 
places, which we will mark *. 

We must now assign B. 1-8^ C, D. 1-2 their actual values, pair the 35 
orders 1,1-6, 11,1-5, 111,1, V,1.2, VI, 1-2, VII, 1-3, VIII, 1-2, IX,1.2, X,l-4, XI, l-^ 
to 18 words, and then assign I,1-XI,8 their actual values. These are expressed in 
this table: 



1,1.6 


- 2' 


X,l-4 


7'- 9 


VIII, 1-2 


16 - 16' 


11,1-5 


3 - 5 


XI, 1-8 


9'- 13 


IX, 1-2 


17 - 17' 


111,1,* 


5'- 6 


V,l-2 


13'- 14 


B. 1-8 


18 - 25 


VI, 1-2 


6'- 7 


VII, 1-3 


14'- 15' 


C 
D.1-2 


26 

27 - 28 



Now we obtain this coded sequence: 






18 , 


27 S 


10 


23 h- 


25 h- 


20 


5o 


1 


19 , 


26 S 


11 


e Cc, 


26 


21 


2/3N 


2 


20 , 


5 Sp 


12 


25 h 


26 S 


22 


13 


3 


26 , 


4 Sp 


13 


3 C 


L 


23 


(p^-l>o 


4 


- R 


, 21 X 


14 


28 S 


26 


24 


6„ 


5 


- C 


R 


15 


23 h- 


, 17 Cc 


25 


U 

lo 


6 


, 


28 S 


16 


24 


7 C 


26 


U 


7 


22 , 


5 Sp 


17 


20 


7C' 


27 


.. 


8 


27 , 


28 h 


18 







28 


.. 


9 


27 S 


, 26 


19 


P« 









Hie' durations may be estimated as follows: 

I: 225 U, II: 270 \i, III: 30 ^i, V: 55 |i, VI: 75 [i: VII: 125 \x, VIII: 75 [i, 
IX: 75 \i, X: 150 p., XI: 300 |jl. 

Total: I + II X (N+l) + III x 2 + VI x (K + i) + (V + VII) x ^ + 

2 2 

+ VIII X (N . 1 ) + IX + (X + XI) x(N+l) = 

= (225 + 270 (N+l) + 60 + 75 (5i! + l) + 180 ^ + 75 (^ - 1) + 75 + 450 (N+D) \x 

2 2 2 

= (885 N + 1080) g.^ (.9 N + 1.1) m. 



10.4 We evaluate next by the formula (9), and this time we keep yJ^, k 
general. We had Z, k = 0, 1, . . . ,N and k-^>0. This excludes k = as well as ^If = N. 
It is somewhat more convenient to write Z-1 f or /f. Theuy?, k = /i*, . . . ,N and k-^^0. 
Thus (9) becomes 



(10) 



k44 p 



N 



f(z)dz = [f(^) + f(^) + ... + f(^) ] i + 

N N N N 

+ [f(^) . £(^) - f(5^) + i(^)] ri- + cp2 , 

N N N N 24N ^^ ' 



j92| i 



Mf" 



180 



k-Jf < Mf"" 1 

N6 180 N'^ 



Finally the requirement k-/^0 may be dropped, since for k-X^ 





Thus N andyf, k - 1, . .. ,N are subject to no further restrictions. 
We state accordingly: 
PROBLEM 11. 
Same as Problem 10, with this change . It is desired to evaluate the integral 

f(z)dz (/, k = 1,...,N) by means of the formula (10) (for k ^ >?, f or k < ^ 

interchange Jf, k). 

The expression to be computed is 



k4^2 

& 

N 



/' 



_ + 



J' for 



k S^ 
k <r 



where J' is defined as follows: Put 



>■■{ 



r = 



for 



k <2' 



then 



and 



24N N N N.N 



J" 



1 k'-l ,u, 
i 2^. f(ii) 

N h=^' N 



The last equation amounts to the inductive definition 



''^'-l 







\ - \.i ^ h f 4^ . 



N N 



J" 



'k'-l 



The storage areas A, B, C, D and the induction index h will be treated the 
same way as in Problem 10, but h runs now over ^\ /If'+l, ..., k^-1. 

We can now draw the flow diagram as shown in Figure 10.2. It will be found 
convenient to store the contents of E. 1-2, (p+k') , (ptJ^O , in the same place where 
(p+k)., (p+j!*) are originally stored, which proves to be B. 1-2. This simplifies II 
and reduces the memory requirements, but since we wish to have B. 1-2 at the end of 
the routine in the same state in which they were at the beginning, it requires restor- 
ing B. 1-2 in DC. VIII, on the other hand, turns out to be entirely unnecessary. 
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£.1 (/\.k% 
2 f/)./^« 



(Tv M?^ 




M. A'A," 01.^1 <o£^ 







C (A.l% 
D O 



M 



«*a* 



to D 



i'-*h 



K+l-^h 



m 



<2 (i4.A«;, 



-fr 






i^S 



"^p^ 
B 



XL 



D /'^^^c. 



h«k' 



6 



■© 






sa 



@^ ^-^''-2^P(4v-^r/;]^^[f(^)-f^t.i>4:0 



^•'^kr^i 



K,i 
^ 



D tT 



m 



@-»— J"»J^'toi> 



s • 



(2^-»^,rs-/'*.j) 




JE 



Figure 10.2 
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The static coding of the boxes I- DC follows: 



B.1 
2 

I.l 
2 
3 

B.3 

11,1 

2 

( 

111,1 

2 

3 

4 

5 

6 

B.4 

111,7 

8 

( 

IV,1 

2 

B.5' 
IV, 3 

4 
( 

V,l 
2 

3 

( 

VI, 1 

2 
3 

[ 

B.6 

VI, 4 

5 

6 

7 

B.7 



(to III.l) 



(to IV, 1) 



(to IV, 1) 



(to V,l) 



(to VI.l) 



(p+k)o 




(P*^>0 




B.1 




B.2 


h- 


II. 1 


Cc 


(ai)o 




B.3 




VII, 23 


Sp 


B.1 




s.l 


s 


B.2 




B.1 


s 


s.l 




B.2 


s 


(012)0 




B.4 




VII,23 


Sp 


B.2 




C 


s 







B.5 




D 


s 


C 




B.1 


h- 


VII, 1 


Cc 


C 




VI, 3 


Sp 


- 


R 


p+h 


R ] 


1 




N 




B.6 


X 


D 


h 


D 


S 


C 




1 





Ac 
Ac 



R 

Ac 

Ac 

D 
Ac 



(p+k)^ 



Ac 


(ai)o 


VII,23 


«! 


Ac 


(p+k)^ 


s.l 


(p^k)^ 


Ac 


(p^^)^ 


B.1 


(p^l)^ 


Ac 


(p.k)^ 


B.2 


(p^k)^ 


Ac 


(a2)o 


VII,23 


tta 


Ac 


(p+r), 


C 


(ptf-) 


Ac 





D 





Ac 


(p+h)^ 


Ac 


(h-k')^ 


Ac 


(p+h) 

^ 


VI,3 


p+h 



R 



f(!i) 

N 



1 f(h) 
N N 

•^h " \-l 



1 f(h) 

N N 



(p+h). 



.12- 



VI, 8 

9 

( 
VII, 1 

2 
3 

4 
5 



[ 

7 

8 

B.8 

VII, 9 

10 

11 

12 

13 

14 

15 
16 

[ 

17 
[ 

18 

19 

20 

21 



22 
23 



(to V,l) 



B.7 


h 


C 


S 


B.2 




VII, 6 


Sp 


B,7 


h- 


VII, 5 


Sp 


p+^'-l 




.. 


h- 


p+r 


h- 


s.l 


s 


s.l 


R 


1 




24N 




B.8 


X 


D 


h 


D 


S 


B.l 




VII, 17 


Sp 


B.7 


h- 


VII, 16 


Sp 


p+k'-l 


- 


-_ 


h 


p+k- 


h 


s.l 


s 


s.l 


R 



B.8 



D 



D 



S 

c 

c ] 



Ac 
C 

Ac 

VII, 6 
Ac 

VII, 5 

Ac 

Ac 

s.l 

R 

Ac 
Ac 
D 

Ac 

VII, 17 

Ac 

VII, 16 

Ac 

Ac 

s.l 

R 

Ac 

Ac 

D 



(p+h+1)^ 

(p+h+1) 

^ o 

p+J" 

(p+ir'-i)^ 
p+^'-l 

N 

f(£lzl) . f(il) 

f(A) . f(ii) 

N N 

£(£1^1) . f(r ) 

N N 



h- 



1 
24N 



iiiiLzl) . i{iL)) 



N 



N 



J" +a_ (f(^:jd) . £(11)) 

24N N N 

J" + _i_ (f(i^) . i{lL)) 

24N N N 

(p+k')^ 

p+k' h 

(p+k'-l) 

p+k'-l 

-f(i^) 

N 

fcki) . f(k:^) 
f(kl) . f(iA) 
f(kl) . i(£i) 

N N 

J- (f(!L!.) . f(kLJ:) 

24N N N 



J" + 



2m 



N 



f(i^)) 

N 



24N N N 
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VIII 




-. 




(to 


e) 






IX,1 




D 


. 


2 




D 


S 


3 




B.l 




4 




s.l 


s 


5 




B.2 




6 




B.1 


s 


7 




s.l 




8 




B.2 


s 


(i 


to e) 







Ac 


J = -J' 


D 


J 


Ac 


(p+i)^ = (p+k') 


s.l 


(p.^)° 


Ac 


(p+k)^= (p+r) 


B.l 


(p+k)^ 


Ac 


(p.^)^ 


B.2 


(p+jr) 



Note, that the box VIII required no coding, hence its inmediate successor 
(e) must follow directly upon its immediate predecessor. However, this box has actually 
no immediate predecessor; VIII, 1 corresponds to a^, which may appear (by substitution at 
11,2) in the C order VII, 23. Hence VIII, 1 must be replaced by e in 0L±.^ 

The ordering of the boxes is I, III, IV, V, VI; II; VII; IX (VIII omitted, 
cf. above), and IV, V, e must also be the immediate successors of II, VI, IX, respectively. 
This necessitates the extra orders: 



11,3 


IV, 1 


C 


VI, 10 


v,i 


C 


IX, 9 


e 


C 



Oi± corresponds to VIII, 1, i.e., to e (cf. above), Oi^ corresponds to IX, 1. 
This implies, as in the corresponding situation in EVoblem 10, that IX, 1 and e must have 
in the final enumeration numbers of the same parity. (3 need not be considered, since it 
represents a fixed remote connection and therefore does not appear in the above detailed 
coding. 

We must now assign B. 1-8, C, D, s.l their actual values, pair the 63 orders 
1,1-3, 11,1-3, 111,1-8, IV, 1-4, V,l-3, VI, 1-10, MI, 1-23, IX, 1-9 to 32 words, and then 
assign I,1-IX,9 their actual values. Ihese are expressed. in this table: 

39 



1,1-3 


- 1 


VI, 1-10 


9 - 13' 


B.1-8 


32 


III. 1-8 


1'- 5 


11,1-3 


14 - 15 


C 


40 


IV, 1-4 


5'- 7 


VII, 1-23 


15'- 26' 


D 


41 


V,l-3 


7'- 8' 


K,l-9 


27 - 31 


s.l 


42 



e is supposed to have the parity of IX, 1, i.e., to be unprimed. 
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Now we obtain this coded sequence: 






32 , 


33 h- 


1 


14 Cc, 


32 


2 


42 S , 


33 


3 


32 S , 


42 


4 


33 S . 


35 


5 


26 Sp', 


33 


6 


40 S , 


36 


7 


41 S . 


40 


8 


32 h-, 


15 Cc' 


9 


40 , 


10 Sp 


10 


- R, 


37 X 


11 


41 h , 


41 S 


12 


40 , 


38 h 


13 


40 S . 


7C' 



14 


34 , 


26 Sp' 


15 


5 C 


33 


16 


18 Sp, 


38 h- 


17 


17 Sp', 


- 


18 


- h- 


42 S 


19 


42 R 


, 39 X 


20 


41 h 


41 S 


21 


32 


23 Sp' 


22 


38 h- 


23 Sp 


23 


- - 


- h 


24 


42 S 


42 R 


25 


39 X 


41 h 


26 


41 S 


- C 


27 


41 - , 


41 S 



28 


32 , 


42 S 


29 


33 , 


32 S 


30 


42 . 


33 S 


31 


eC , 


_. 


32 
33 
34 


(p+k)o 

^0 




35 


27o 




36 
37 

38 
39 



1 

N 
24N 




40 


-- 




41 






42 


-. 





Tlie durations may be estimated as follows: 

I: 125 \i, II: 125 jj., Ill: 300 g., IV: 150 |i, V: 125 \x, VI: 445 |i, VII: 1015 |i, IX: 350 ^l. 

Total: I + (II or III) + IV + V x (|k-l| + 1) + VI x | k-Jf| + VII + (VIII or IX) 

maximum = (125 + 300 + 150 + 125 (\k-Z\ +'-!) + 445 | k-i\ + 1,015 + 350 |i = 

= (570 \k-y!C\ + 2,065) n 
maximum = (570 (N-1) + 2,065) \i = (570 N + 1,495) usi (.6 N + 1.5) m. 

10.5 We pass now to the interpolation problem. 

Lagrange's interpolation formula expresses the unique polynomial P(x) of 
degree M-1, which assumes M given values pj^,. ..,p[^^ at M given places x]^, • . . ,Xj^j, 
respectively: 



(I) 



M 



P(x) - P(xi, pi;...; 5^, Pm|x) = .1^ Pi 



-S^i^i) (x-x,) 



izl 



Mr..i) 






(Xi-Xj) 



There would be no difficulty in devising a (multiply inductive) routine 
which evaluates the right hand side of (1) directly. This, however, seems inad- 
visable, except for relatively small values of M. Hie reason is that the denominators 



?*/ are li 



j=l 



are likely to 'be inadmissibly small. 
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This point may deserve a somewhat more detailed analysis. 

From our general conditions of storage and the speed of our arithmetical 
organs one will be inclined to conclude that the space allotted to the storage of the 
functions which are evaluated by interpolation should (in a given problem) be comparable 
to the space occupied by the interpolation routine itself. TTie latter amounts to about 
100 words. (Problem 12 occupies together with Problem 13. a or 13. b or 13. c 99 or 101 
or 106 words, respectively, cf. 10.7 or 10.8 or 10.9, respectively. Other possible 
variants occupy very comparable amounts of space.) One problem will frequently use 
interpolation on several functions. It seems therefore reasonable to expect that each 
of these functions will be given at ^ i • 100^^25 points. (lliis means that the N of 
10.7 will be^^ 2^ -- not our present M ! Note also, that the storage recifuired in this 
connection is N in Problem 13. a, but 2N in Problems 13. b and 13. c.) Hence we may 
expect that the points X]^,...,xiL« will be at distances of the order '-^ 2'^ between 
neighbors. 



Hence I ^i-Xj | ^ 



^jM 



i=l 



1-J 



(X£-X-) 



and 



- (i-l) ! (M-i) I 2-5^^-1^ 



The round-off errors of our multiplication introduce into all these products absolute 
errors of the order 2"^^° . Hence the denominators -^-j^ (x^-x-), and with them the 



corresponding terms of the sum Z. in (1), are affected with relative errors of the 

i-l 

order 2-^0; (i-l) ! (M-i) ! 2-^^^-^^ - ' ^ (^'i). The average affect of these 

relative errors is best estimated as the relative error 




2 5M-45 



if.\n 



1^2 



(m-2\ 



,5h4-45 



(M-D! M 




,ai-45 



(M-1) ! vW 



(2(M-1))! 



[(M-1)!? 



(2(M-1))! 



M 



[(M-1)!P 



,5M-45 



On the other hand, an interpolation of degree M-1, with an interval 
length -^^ 2'^ , is likely to have a relative precision of the order ~' C'2" , where 
C is a moderately large number. (The function that is being interpolated is assumed 
to be reasonably smooth. ) 
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Gxisequently the optimum relative precision e which can be obtained with 
this procedure, and the optimum M that goes with it, are determined by these conditions: 



-/ 



(2<M - 1))! 



From this 



%=p¥ 



m 



[(M-l)!] 



1 25M-45 ^ c • 2-^ 



[(M-l)!? 



2lOM-45^ C. 



Now we have 



M 


3 


4 


5 


6 


7 


% 


8-10-« 


10-2 


5 


1.8-103 


6-10^ 



The plausible values of C are in the neighborhood of M = 5, while M = 6 is somewhat 
high and M = 4 and 7 are extremely low and high, respectively. Hence under these 
conditions M = 5 (biquadratic interpolation) would seem to be normally optimal, with 
M = 6 a somewhat less probable possibility. We have 



M 



L4-10-^ 



1.7-10-6 



It follows, therefore, that we may expect to obtain by a reasonable 
application of this method relative precisions of the order --^ 10"® ^^2"^°. This is, 
however, only the relative precision as delimited by one particular source of errors: 
The arithmetical (round-off) errors of an interpolation. Ihe ultimate level of pre- 
cision of the entire problem, in which this interpolation occurs, is thei'efore likely 
to be a good deal less favorable. 

These things being understood, it seems likely that the resulting level 
of precision will be acceptable in many classes of problems, especially among those 
which originate in physics. There are, on the other hand, numerous and important 
problems where this is not desirable or acceptable, especially since it represents 
the loss of half the intrinsic precision of the 40 (binary) digit machine. It is 
therefore worthwhile to look for alternative procedures. 

The obvious method to avoid the loss of (relative) precision caused in 

M(j/i) 

the formula (1) by the smallness of a denominator il (x'-x-), is to divide 

j=l ^ J 

by its factors x--x- (j = 1,...,M and j ^ i) singly and successively. This must. 
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however, be combined with a similar treatment of the numerators TT (x-x-), since 

j=l J 

they may cause a comparable loss of (relative) precision by the same mechanism. A 
possible alternative to (1), which eliminates both sources of error in the sense indi- 
cated, is 

(2) P(x) = .2, Pi TT^ ^^ 



1 



=1 '-I j=l ^i--j 



(2) involves considerably more divisions than (1), but this need not be the dominant 
consideration. Tliere exists, however, a third procedure, which has all the advantages 
of (2), and is somewhat more easily handled. Besides, its storage and induction prob- 
lems are more instructive than those of (1) or (2), and for these reasons we propose 
to use this third procedure as the basis of our discussion. 

This procedure is based on A. C. Aitken's identity 

(3) P(xi, pi;,..; xj^, j^ |x) = -|-^ ^^^2' P2'--- ^' RW 1^^ ■" 



X -X| 



P(xi, Pi;...; x|^.i, v^^i Ix) 



Since (3) replaces an M-point interpolation by two M-1 point interpolations, it is 
clearly a possible basis for an inductive procedure. It might seem, however, that 
the reduction from an M-point interpolation to one-point ones (i.e. to constants) 
will involve 2^"^ steps, which would be excessive, since (2) is clearly an M(M-l) 
step process. However, (3) removes either extremity (xj or xj^) of the point system 
X]^, — ,xj^; hence iterating it can only lead to point systems Xj^,...,x- (i, j = 1, . . . ,M, 

i ^ j, we will write j = i+h-1), of which there are only Mi^^ ; i.e., Mifikil , not 

counting the one-point systems. Hence (3) is likely to lead to something like em 

MiMzi) step process. 

Regarding the sizes we assume that xi,...,xj^ as well as x and Pl,---iPM 
lie in the interval -1, 1. We need, furthermore, that the differences X£4.j^-X£, 
X-X-, Pi+h'Pi ^^^^ ^^® ^^ ^^^ interval -1, 1, and it is even necessary in view of 
the particular algebraical routine that we use, to have all differences Pi+h"Pi 
(absolutely) smaller than the corresponding x^+j^-x^. (Cf. VII, 19. ) I.e., we 
must use appropriate size adjusting factors for the p^ to secure this "Lipschitz 
condition". The same must be postulated for the differences P£+i(x)-P^(x) of the 
intermediate interpolation polynomials. All of this might be circumvented to 
various extents in various ways, but this would lead us deeply into the problems 
of polynomial interpolation which are not our concern here. 

We now state: 
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PROBLEM 12. , 

The variable values X]^, ,xj^| (x|< xo"^ ...^ x^) ^^^ the function values 

P1»'**»PM ^^^ stored at two systems of M consecutive metnory locations each: 
q, q+1, ..., q-*-M-l and p, p+1, . .. , p+M-1. The constants of the problem are p, q, 
M, X, and they are stored at four given memory locations . It is desired to inter- 
polate this function for the value x of the variable, using Lagrange's formula. 
The process of reduction based on the identity (3) is to be used. 

It is clear from the previous discussion, that we have to form the 
family of interpolants 



(4) P.(x) =P(x., Pi;...; x.+j^,i, Pi+h-l 1^) for i = 1,...,M, 

h = l,...,M-i+l . 

Applying (3) topCx^, p^;...; X£4h' Pi+h 1^^ instead of Plx^, Pi;...; xj^, Pm 1^) 
gives 

h+1, V X -X- h , , X;44,-X h, . 

1 ^i-hh-^i ^^^ ^i+h^i ^ 

or equivalently 



Combining this with 

(6) p|(x) = p- 

anct 

M 

(7) P(x) = P(xi, px;...; xjyj, pj^i !x) = Pj^x) 

we have a (doubly) inductive scheme to calculate P(x). 

, Let A and B be the storage areas corresponding to the two intervals of 
locations from p top-«-M-l and from q toq-^-M-l. In this way their positions will be 
A.1,...,M and B. 1,...,M, where A.i and B. i correspond to p+i-1 and q+i-1, and store 
p- and X-, respectively. As in the Prohibits 3, 10, and 11, the positions of A and B 
need not be shown in the final enumeration of the coded sequence. 

The primary induction will clearly begin with the Pii..«fP|^j, i.e., 
Pj(x), . . . ,'P|^(x), stored at A, and obtain from these Pj(x), . . . ,F^_^x(x); then from 
these the P-^ix) , . . , ,P^^2^x) ; from these the P^lx), . . . ,P|^^3(x); etc., etc., to 
conclude with PjCx), i.e., iwith the desired P(x). Tliese successive stages corres- 
pond to h = 1, 2, 3,. . . ,M, respectively. This h is the primary induction index. 
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, - In passing from h to h+1, i.e., from Pi(x), ..., %.k+]^x) to 
P- (x),..., Pjyj^j^(x), the P- (x) have to be formed successively for i ~l,...,M-h. 

This i is the secondary induction index. 

h h . . 
Ihere is clearly need for a storage area C which holds P^^x), . . . ,P|^^j^^j(x) 

at each stage h. As the transition to stage h+1 takes place, each P^(x) will be re- 
placed byl^ (x), successively for all i = 1, ..., M-h. Hence the capacity required 

for C is M-h+1 during the stage h, but for h = 1 the Pi(x) = p^ are still in A, and 
need not appear in C. Hence the maximum capacity for C is M«-l. Accordingly, let C 
correspond to the interval of locations from r to r"HV1-2. In this way its positions 
will be C. 1,...,M-1, ivhere C.i corresponds to r+i-1, and stores P^^x). The positions 
of the area C need not be shown in the final enumeration of the coded sequence. All 
that is necessary is that they be available (i £ . empty or irrelevantly occupied) when 
the instructions of the coded sequence are being carried out by the machine. 

As stated above, the P^(x) occupying C.i have necessarily h ^ 1. To give 
more detail: Pj^(x) moves into the location C.i at the end of the i-th step of stage 
h-1, and remains there during the remainder (the M-h+l-i last steps) of stage h-1 and 
during the i first steps of stage h. 

Further storage capacities are required as follovs: The given data (the 
constants) of the problem, p, q, r, M, x, will be stored in the storage area D. (It 
will be convenient to store them as p^, (q-p-l)^, r^, (M-D^, x. ) Storage will also 
have to be provided for various other fixed quantities (1^, the exit- locations of the 
variable remote connections), these too will be accomnodated in D. Next, the two in- 
duction indices h, i, will have to be stored. As before, they are both relevant as 
position marks, and it will prove convenient to store in their place (M-h)^, (p+i)^ 
(i.e., (A.i+1)^). These will be stored in the storage area E. Finally, the quanti- 
ties which are processed during each inductive step will be stored in the storage area 
F. 

We can now draw the flav diagram, as shown in Figure 10.3. Tlie variable 
remote connections OL and 3 s^re necessary, in order to make the differing treatments 
required for h = 1 and for h > 1 possible: In the first case the P^(x) are equal to 
pj and come from A, in the second case they come from C (cf. above). It should be 
noted that in this situation our flow diagram rules impose a rather detailed showing 
of the contents of the storage area C. 

The actual coding ccntains a number of minor deviations from the flow 
diagram,) inasmuch as it is convenient to move a few operations from the box in which 
they are shovn to an earlier box. Since we had instances of this already in earlier 
problems, we will not discuss it here in detail . On the other hand, sane further 
condensations of the coding, which are possible, but deviate still further from the 
flow diagram, will not be considered here. 
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©^ 



(M-1)^ to £4 



£.1 (M-l)^ 



£A (M-hl 



C.l 



n <^>\ CLYfd h>i 



E.l (M-hh 



-^ 



1 -^A 



3 



6^.2^*0 £.£ 



^-^ 



© 



2 



© 



C.i'' 






•A+i 



Oihe 



E.i (M-hh 
2 (A.Ut)^ 



^^ ^ — ^^ ^^y rTTTTd £• 

©tprrrw ^^rx; «o ^.i '^ , — — _, — . — , 
4 1 j^ l ^Li __ 1 E.l (M^h)^ F.I P^CX) 

^^-^ ' i - ' I i i " ! .i« i ii " i > ...f ir I'liii •mill Ml n i 9 X.. ' f f I 



F.I Ff(X) 
2 — 



®iilp 



M 



J^^ 



- — 1 > — " 



"HT 



" JH 



T(*h'Tt(^*^[^yn(4 



to a 



Ci'pfidf) for i'i M-h 






pf^'M foi^l'<i 






k-^l 



iz 



TTyir 



cr 









1 



Oih€ 



(i-fh-M)^ 



(M-h-jt). 



14 



+'\ 



I 
I 



tz 



# 



hf 



i ^ M-h 



/4 



/;3 



TZ 



E.l(M-h)o 
Z(A.i*Z\ 


F.I Pi'^^ (X) 


fU ■ 



T\^^ (K) to FA 






1-9-1-^1 



8 



O 






/«' 



I 
I 



A + i-^ h 



<D 



E.l(M-k-ll 



# 



A,i 



/« 



CI P(X}'PJX) 



A,i 
B.i 



<i) 






E.l rM-h)e 

C.i Pi(x) fofltM-hfl 



IG-URE 



10.3 
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Tlie static coding of the boxes I-XI follows: 



D.l 


(M.l)„ 




I.l 


D.l 




2 


E.1 


S 


D.2 
1,3 
4 
D.3 
1,5 


(ai)o 
D.2 
11,4 
Oi)o 
D.3 


S] 


6 
(to II, 1) 
D.4 


111,6 

Po 


Si 


5 
II.l 


D.4 




2 


D.5 


h 


3 


E.2 


s 


4 


.. 


c 


111,1 


D.4 


c 


2 
3 


111,3 

P 


S] 


4 


F.l 


s 



D.6 

IV, 1 

2 
3 



V.l 
2 



(to 111,4) 



E.2 



D.6 
IV, 3 



V,2 



VI. 1 
2 

3 

4 



(to VII, 1) 



r+i 
F.2 



C 

c ] 



Sp 



Sp 



p+i 




F.2 


S 


D.4 


h- 


D.6 


h 


VI, 4 


Sp 



Ac 
E.l 

Ac 
11,4 

Ac 
111,6 



Ac 



Ac 
IV, 3 

Ac 

V,2 

Ac 

F.2 

Ac 
Ac 
VI, 4 

Ac 

F.2 



(M-l)j 
(M-1), 



(ai)( 



(3i)o 
3. 



Ac 


Po 


Ac 


(p-^Do 


E.2 


(p+l)o 


Ac 


Po 


111,3 


p 


Ac 


Pj(x) = pj 


F.l 


Pj if reached from 111,3 



P|(x) if reached from IV 

(p+l)_if reached from 111,4 
(p+i)Qif reached from IX (i.e. 



p5(x) 

p+i 

Pi-Kl 
Pi+1 



"0 

(r+i)^ 
r+i 



F*^l(x) 
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VII.l 
2 
3 
4 
5 
6 
D.7 
VII, 7 
8 
9 

10 
11 
12 
13 
14 
[ 

15 

[ 

16 

17 

18 

19 

D.8 

VII, 20 

21 

[ 



(to VII, 1) 



22 
23 

24 

25 



E.2 




D.4 


h- 


D.6 


h 


D.5 


h- 


VII, 25 


Sp 


E.2 




(q-p-l)o 




D.7 


h 


VII, 15 


Sp 


VII, 21 


Sp 


D.l 


h 


E.l 


h- 


D.5 


h 


VII, 14 


Sp 



q+i+h-1 

q+i-1 
s.l 

F.2 

F.l 



r+i-1 



h- 
h- 

S 

h- 



X 






D.8 






— 


h- 




q+i-1 


h- 


] 


s.l 


S 




s.l 


X 




F.l 


h 





s 
s 



(to VIll.l) 



Ac 


(p^)o 


Ac 


io 


Ac 


(r+i). 


Ac 


(r+i-l)„ 


VII, 25 


r+i-l 


Ac 


(P^i)o 


Ac 


(<l+i-l)o 


VII, 15 


q+i-1 


VII,21 


q+i-1 


Ac 


(q+i4M-2)^ 


Ac 


(q+i+h-2)^ 


Ac 


(q+i+h-Dg 


VII, 14 


q+i+h-1 



Ac 

Ac 
s.l 

Ac 

Ac 

R 

Ac 

Ac 

s.l 

Ac 
Ac 



C.i 



4+h 



''i+h-'^i 
''i+h-'^i 
pS,i(x) 



F*+l('')-Pi<*> 
P^.n(x)-P^(x ) 



X-Xl 



x-x. 



1 

X -X 



h- 
h- 



i (I*+i(x)-P^(x)) 



"i+h-'^i '"i+l'"' -1^ 
I**^(x) = P^(x) + 



*i+h "i ^^ ^ 



pJ^Mx) 
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VIII, 1 




E.2 




2 




D.5 


h 


3 




E.2 


s 


4 




■ E.l 


h- 


,5 




as 


h- 


6 




D.4 


h- 


7 
(to 


IX, 1) 


X.l 


Cc 


IX. 1 




F.2 




2 
(to 


111,5) 


F.l 


s 


X.l 




E.1 




2 




as 


h- 


3 




E.l 


S 


4 




as 


h- 


5 
(to 
D.9 


e) 


XI, 1 

(«2)0 


Cc 


XI, 1 


V 


a 9 




2 

a 10 




11.4 

(|32)0 


Sp 


XI, 3 




a 10 




4 
(to 


II, 1) 


111,6 


Sp 



Ac 


(p+i)o 


Ac 


(p-i-1). 


E.2 


(pH+l)„ 


Ac 


(p+i+h-M+l)g 


Ac 


(p+i-Hi-M) ^ 


Ac 


(i+h-M)„ 


Ac 


p5,i(x) 


F.l 


p5,i(x) 


Ac 


(M.h)o 


Ac 


(M-h-1)^ 


E.l 


(M-h-l)„ 


Ac 


(M-h-2)„ 



Ac 


(«2)o 


II.4 


02 


Ac 


(P2)0 


111,6 


fe 



The ordering of the boxes is I, II; III; IV; V, VII, VIII, IX; X; XI; 
VI, and VII, e, II must also be the immediate successors of VI, X, XI, respectively, 
and 111,4 and 111,5 must be the immediate successors of IV and IX, respectively. 
This necessitates the extra orders 



VI, 6 


VII, 1 


C 


X,6 


e 


G 


XI, 5 


II. 1 


C 



and 



IV, 4 
IX. 3 



111,4 
III. 5 



C 
C 



^it °ht Pi. Pa correspond to 111,1, IV, 1, V,l, VI, 1, respectively. 
Hence in the final enumeration 111,1, IV, 1 must have the same parity, and V, 1, 
VI, 1 must have the same parity. 
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We must now assign D.l-lO, E.1-2, F. 1-2, s.l their actual values, pair 
the 75 orders 1,1-6, 11,1-4, 111,1-6, IV, 1-4, V,1.3, VI, 1-6, VII, 1-25, VIII, 1-7, 
IX, 1-3, X, 1-6, XI, 1-5 to 38 words, and then assign I,1-XI,5 their final values. 
Ihese are expressed in this table: 



I, 


1-6 




-2' 


VII, 


1-25 


11 '-23' 


VI, 1-6 


35 -37' 


11. 


1-4 




3 -4' 


VIII, 


1-7 




24 -27 


D.1-10 


38 -47 


in, 


1-6 




5 .7' 


IX,l-3 




27 '-28' 


E.1-2 


48 -49 


IV, 1-4 




8 -9' 


X,l-6 




29 -31' 


F. 1-2 


50 -51 


V,1.3 




10 -11 


XI. 


1-5 


* 
> 


32 -34' 


s.l 


52 






Now we obtain this coded 


sequence: 











38 , 


48 S 


18 






, - h- 


36 


36 Sp', - 


1 




39 , 


4 Spt 


19 




52 S , 


51 


37 


51 S , 11 C 


2 




40 , 


7 Sp' 


20 




50 h-, 


52 ^ 


38 


(M.l)„ 


3 




41 , 


42 h 


21 




45 , 


- h- 


39 


5„ ° 


4 




49 S , 


- C 


22 




52 S , 


52 X 


40 


A 


5 




41 ■ 


6 Sp 


23 




50 h 


- S 


41 


Po 


6 




, 


50 S 


24 




49 , 


42 h 


42 


1° 


7 




49 , 


- C 


25 




49 S , 


48 h- 


43 


u 
r 


8 




43 . 


9 Sp 


26 




42 h-, 


41 h- 


44 


(q-p-l)o 


9 




, 


6C' 


27 




29 Cc 


51 


45 


X 


10 




10 Sp' 


, - 


28 




50 S , 


7 C 


46 


8_ 


11 




51 S , 


49 


29 




48 


, 42 h- 


47 


A 


12 




41 h-, 


43 h 


30 




48 S 


42 h- 


48 


u 


13 




42 h-, 


23 Sp' 


31 




32 Cc 


, eC 


49 


- 


14 




49 , 


44h 


32 




46 


4Sp' 


50 


. 


15 




18 Sp', 


21 Sp' 


33 




47 


7 Sp' 


51 


- 


16 




38 h , 


48 h- 


34 




3 C 


- 


52 


- 


17 




42 h , 


18 Sp 


35 




41 h- 


, 43 h 







The durations may be estimated as follows: 

I: 225 \x, II: 150 n, III: 225 \i, IV: 275 ji, V: 125 \i, VI: 225 \i., VII: 1140 ^, 
VIII: 275 \i, JX: 200 ^, X: 225 p., XI: 200 p. 

Total: I + II X (M-1) + III + IV x (M.2) + V x (M-1) + VI x (M-lKM-2) + 

+ (VII + VIII) M(M-l) + IX X (M-l)(M-2) + X X (M-1) + XI x (M-2) = 
2 2 

= (225 + 150 (M-1) + 225 + 275 (M-2) + 125 (M-1) + 225 (M-lMM-2) + 

+ 1.415 M(M-1) + 200 (M.l)(M-2 ) + 225 (M-1) + 200 (M-2)) [i = 
2 2 



= (920 fcP - 370 M - 575) \i »(.9M2 - .4 M - .6) m. 
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10.6 We now pass to the problem of interpolating a tabulated function 
based on the coding of Lagrange's interpolation formula in Problem 12. 

PROBLEM 13. 

The variable values y-j^, ...,yjy^ (y|^< y2'^ .••< yjyj) and the function values 
qi>« -fC^ are stored at two systems of N consecutive memory locations each: 
q. q+1, ... , q+N-1 and p, p+1, . . . ,p4N-l. The constants of the problem are p, q, N, M, 
y, to be stored at five given memory locations. It is desired to interpolate this 
function for the value y of the variable, using Lagrange's formula for the M points 
y- nearest to y. 

Hie problem should be treated differently, according to whether the 
yi#»«-»ytS[ are or are not to be equidistant. 

PROBLEM 13. a. 

yi,...,>^ are equidistant, i.e., y^ = a + ^^^^ (bra). In this case 
only y]^ = a and Ypj ~ b need be stored. 

PROBLEM 13.b. 

yi»*'*>7N» ^^^ unrestricted. 

In both cases ou!r purpose is to reduce FVoblem 13 to Problem 12, with 

xj,...,xj^l equal to 7^* •• • »yk-+M-l' ^^ Pl» • • • »Pm ^^^^^ ^^ %' • ' * »^k+M-l' ^^^^^ 

k = ] , . . . ,N-M+1 is so chosen that the yui.. •tyk'+M-l ^^^ ^^ close to y as possible. 

10.7 We consider first Problem 13. a. 

In this case the definition of k, as formulated at the end of 10.6, 
amounts to this: The remoter one of yi = a "^ ^^^— (b-a) and yujju i = a + '^ — ^ (b-a) 

should lie as close as possible to y. This is equivalent to requiring that their 
mean, 9(yk"^k-%l-l^ "= a + -^ ■ (b-a), lie as close as possible to y, i.e.,. that 

k lie as close as possible to (N-1) ^^^ - '^^ , 

b-a 2 

There are various ways to find this k = R, based on iterative trial and 
error procedures . Since we will have to use a method of this type in connection with 
Problem 13. b, we prefer a different one at this occasion. 

lliis method is bas<id on the function Iz), which denotes the integer 
closest to z. Putting 



k* = /(N-1) i^ - Mr3l ^ 
1^ b-a 2 J 



we have clearly 
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{ 



= 1 for k* g 1, 

= k* for 1 ^ k g N-M+1, 

= N-M+1 for k* ^ N-M+1 



Now the multiplication order (11, Table II), permits us to obtain iz] 
directly. Indeed, the round off rule (cf. the discussion of the order in question) 
has the effect that when a product uv is formed, the accumulator will contain 

uv = 2"^® {2^^uv), and the arithmetical register will contain uv = 2''^uv - {2^,^uv). 

Hence putting u = 2-='«(N-l) and v = i^ - -M=^ will produce "^ ^ 2'^^ {2=^^} = 2-^V 

in the accumulator. 

After k* and E have been obtained, we can utilize the routine of Problem 12 
to complete the task. This means, that we propose to use the coded sequence 0-52 of 
10.5, and that we will adjust the coded sequence that will be formed here, so that it 
can be used in ccmjunction with that one of 10.5. 

Among the constants of Problem 12 only p, q, and x need be given values 
which correspond to the new situation, x is clearly our present y. p,...,p-*M-l are 
the positions of the p^,. . mP[^ o£ Problem 12, i.e., the positions of our present 
^E> • • 'j^g-fM-l* i»e. , they are p-+ic-l, . . . ,p"^E'^-2. Hence p = p"^-l. q, .•.,q+M-l are 
the positions of the x],...,xj^ of Problem 12, i.e., the positions of our present 
y^, . . . ,yT^4j^_]^. Since we determined in the formulation of Problem 13. a, that only 
yj = a and )^ - b will be stored, but not the entire sequence yi,-»«>y][^> this means 
that the desired sub-sequence Yr* * • • ^Yj^-fM-l does not exist anywhere ab initio. 

Consequently q may have any value, all that is needed is that the posi- 
tions q,...,q"'^-l should be available and empty (or irrelevantly occupied) when the 
Qoded sequence that we are going to formulate begins to operate. This sequence must 
then form 



and place it into the position q+i-1 for all i - 1,...,M. 

It might seem wasteful to form x- first, then store it at q+i-1, and 
finally obtain it from there by transfer when it is needed, i .e. during the period 
VII, 1-19 of the coded sequence of 10.5. Ohe might think that it is simpler to form 
X£ when it is actually needed, i.e., in VII, 1-19 as stated above; the quantities 

. X -X- 

needed are more specifically x • , x-4j., in the combination -^ ), and thus avoid 

1 in ^i+h"^i 

the transfers and the storage. It is easy to see, however, that the saving thus 

effected is altogether negligible, essentially because the size of our problem is 

proportional to M^ (cf. the end of 10^5), while the number of steps required in 

forming and transferring the x^'s in the first mentioned way is only proportional 

to M. (Remember that M is likely to be ^7, cf. the beginning of 10.5.) It does 

therefore hardly seem worthwhile to undertake those changes of the coded sequence 

of 10.5 which the second procedure would necessitate, and we will adhere to the 

first procedure. 
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In assigning letters to the various storage areas to be used, it must be 
remembered that the coded sequence that we are now developing is to be used in con- 
junction with (i.e., as a supplement to) the coded sequence of 10.5. The latter 
requires storage areas of various types: D-F, which are incorporated into its final 
enumeration (they are 38-51 of the 0-52 of 10.5); A (i.e., p, . . . ,p-^-l), which will 
be part of our present A (ire., of p, . . . ,p-HN-l, it will be p+K-l, . . . ,p"^lc+M-2 ; B 
(i.e., q, . . . ,q-^-l), which may be at any available place; and C (i.e., r, . ,. ,r+M-2), 
which, too, may be at any available place. Therefore we can, in assigning letters 
to the various storage areas of the present coding, disregard those of 10.5, with 
the exception of B, C. Furthermore, there will be no need to refer here to C (of 
10,5), since the coded sequence of 10.5 assumes the area C to be irrelevantly occupied. 
It will be necessary, hov/ever, to refer to B (of 10.5), since it is supposed to con- 
tain the X£ = Yj^+j^.i (i ~ 1,...,M) of the problem. We will therefore think of the 
letters which are meant to designate storage areas of the coded sequence of 10.5 as 
being primed. This is, according to the above, an immediate necessity for B. We 
can now assign freely unprimed letters to the storage areas of the present coding. 

L^t A be the storage area corresponding to the interval from p to 
p-^f^-1. In this way its positions will be A.l,... .,N, where A. i corresponds to 
p+i-1 and stores q^^. As in previous problems, the positions of A need not be 
shown in the final enumeration of the ceded sequence. 

The given data of the problem are p, M, N, a, b, y, also the q, r of 
the coded sequence of 10.5. M, r are already stored there (as (M-D^, r^ at 38, 4S; 
and y, too (it coincides with x at 45). q, however, occurs only in combination with 
Ic (as (q-p-l)^ = (q-p-lc)^ at 44) and p, of course, contains k (as p^ = (p+Ic-D^ at 
41); and E originates in the machine. Hence, 41, 44 must be left empty (or, rather, 
irrelevantly occupied) when our present coded sequence begins to operate, and they 
must be appropriately substituted by its operations, q, however, must be stored as 
one of the constants of our present problem. Hence the constants requiring storage 
now (apart from M, r, y which we stored in 0-52, cf. abcvre) are p, q, N, a, b. 
They will be stored in the storage area B. (It will be ccnvenient to store them as 
a, b, (N-1)q, (p-Dq, (q-l)o- ^^ will also need 2"^® and 1^; we will store the 
former in B and get the latter from 42. ) The induction index i is a position mark; 
it will be stored as (q+i-D^ (i.e. , (B'. i)^) in the storage area C. The quantities 
which are processed during each inductive step will also be stored in the storage 
area C. 

Our task is to calculate K; to substitute p^ = (p+k-D^ and 
(q-p-D^ = (q-p-k)^ into 41 and 44; and then to transfer x- = YF+i-l from A.lc+i-1 
to B'.i. This latter operation is inductive. Finally, the control has to be sent, 
not to e, but to the beginning of 0-52, i.e., to 0. 

We can now draw the flow diagram, as shown in Figure 10.4. 
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5E 
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C.X Xj^i 
2 (B/i^l)o 
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The actual coding will again deviate in some minor respects from the flow 
diagram, as in previous instances. In connection with this we will also need some 
extra storage capacity in C (C. 1.1). 

A matter of somewhat greater importance is this: We have noted before 
that since our machine recognizes numbers between -1 and 1 only, integers t must be 
stored in sane other form. Frequently the form of a position mark, t^t is per se 
more natural than any other (cf . 8.2); occasionally 2'^®Jl can be fitted more easily 
to the algebraical use to be made of t (cf. I in the present coding); sometimes y 
is most convenient. (We could, of course, add to this list, but the three above forms 
seem to be the basic ones .) The transitions between these three forms are easily 
effected by using multiplications and divisions, but it seems natural to want to 
achieve the transitions between the two first ones {t^ and 2"'^®X) in a more direct 
way. 

This can be achieved by means of the partial substitution orders 18, 19 
of Table II (x Sp, x Sp'), if they are modified as indicated in the Remark imne- 
diately preceding this chapter. We propose to use these orders now for arithmetical 
rather than logical (substitution) purposes, for positions x which contain no orders 
at all, but which are storing numbers in transit. Specifically: With t^ in the 
accumulator, x Sp' produces 2'^^t at the position x; with 2"^®X in the accumulator 
X Sp produces 2"^®X at the position x, and a subsequent x h produces / in the 
accumulator. 

We mention that ~ can be obtained from j^ or 2" ^^X by dividing them 
into 1q or 2"^®, respectively! and this without more than the usual loss in precision: 
Indeed, our division p : a is precise within an error 2'^®, no matter how small o is 
(subject, of course, to the condition | a| ^ |p | ), provided that p, are given exactly. 
(If they are not given exactly, then their errors are amplified by ^ , ^, respectively. 
In this case a small O is dangerous, even though | a| > | pi .} In the present case 
pj a are given exactly. Forming 1^, t^t as well as 2"^®, 2"^ I, involves no round- 
offs. 

These methods will be used in our present coding: For transitions 
between t^ and 2-30j:,cf. 1,22-23, 111,3-4; VII, 1-3; For forming reciprocals: -i- 

from (N-1)q as 1^: (N-1) in VIII, 4-5, and the very similar case of ^"^ from 
(M-3)^, (N-Dq as (M-3)j^: 2(N-1)^ in 1,9-15. Regarding -i- we also note this: 

We need ^ in VIII. We will form -i- = 1„: (N-l)„ first and ^ = -i- X (b-a) 

N-1 N-1 ° ° N-1 N-1 

afterwards. Forming 1. x (b-a) first and SiS » [l X (b-a) ]: (N-l)„ afterwards 

would lead to a serious loss of precision, since 1^ X (b-a) plays the role of p 
above, and as it involves a round-off it is not given exactly, and hence may cause 
a loss of precision as indicated there. 

We will have to refer in our present coding repeatedly to 0-52 in 
10.5. It should therefore be remembered that this coded sequence is now supposed 
to be changed insofar that 41, 44 are irrelevant and 45 contains v. 38, 43 contain 
(M-1)q, r^, as in 10.5. 
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The static coding, of the boxes I-X follows: 



B.1 


a 








2 


b 








1,1 


B.2 




Ac 


b 


2 


B.1 


h- 


Ac 


b-a 


3 


s.l 


S 


s.l 


b-a 


4 


45 




Ac 


y 


5 


B.1 


h- 


Ac 


y-a 


6 


s.l 


+ 


R 


yjr.a 

b-a 


7 




A 


Ac 


Xza 

b-a 


8 


s.l 


S 


s.l 


b-a 


B.3 


(N-I), 








1.9 


B.3 




Ac 


(N-l)„ 


10 


B.3 


h 


Ac 


2(N-1), 


11 


8.2 


s 


s.2 


2(N.l) 


12 


38 




Ac 


(M-1), 


13 


42 


h- 


Ac 


(M-2)^ 


14 


42 


fa- 


Ac 


(M-3)„ 



15 

16 
17 

18 

19 

20 
21 
22 
23 

24 

25 



(to 11,1) 



B.4 
11,1 

2 

3 



5.2 



A 

S 



3.1 



h- 



s.l 


S 


s.l 


R 


B.3 




s.2 


Sp 


s.2 


X 


C.l 


S 


2-39 




C.l 




B.4 


h- 


111,1 


Cc 



(to IV, 1) 



R 

Ac 

s.2 

Ac 

Ac 

s.l 
R 
Ac 
s.2 

Ac 
CI 



Ac 
Ac 



M-3 

2(N-1) 

M-3 
2(N-1) 

M-3 

2(N-1) 



(M-3)^ 
2(N.l), 



b-a 




V 


b-a 


. M-3 
2(N-1) 


V 






V 






(N 


-l)o 




u 


= 2-3E 


*(N-1) 


2- 


^\* -- 


'■ uv 


2- 


3^. 




2- 
2- 


39j^* 


-1) 
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111,1 


B.3 




Ac 


(N-l)o 


2 


38 


h- 


Ac 


(N-M), 


3 


42 


h 


Ac 


(N-M+1)^ 


4 


C.1.1 


Sp' 


C.1.1 


2-3 8(N-M+l) 


5 


CI 




Ac 


2-^\* 


6 


C.1.1 


h- 


Ac 


2-3e(k*-N-tM-l) 


7 


v,i 


Cc 







IV. 1 

( 

v.i 

2 
VI, 



(to VI.l) 
L 
(to V,2) 



(to VII, 1) 
(to VII, 1) 



B.4 

C.1.1 
CI 



VII,1 




CI 




2 




s.l 


Sp 


3 
B.5 




s.l 
(P-I)o 


h 


VII, 4 




B.5 


h 


5 




41 


S 


6 




s.l 


S 


B.6 




<q-l)o 




VII,7 




B.6 




8 
9 
(to 
VIII, 1 


VIII, 1) 


s.l 
44 

B.2 


h- 

s 


2 




B.1 


h- 


3 




s.l 


s 


4 




42 





B.3 



6 


s.l 


X 


7 


C.1.1 


S 


8 


C.1.1 


R 


9 


C.1 




10 


B.4 


h- 



Ac 

Ac 
CI 



Ac 

s.l 

Ac 

Ac 
41 
s.l 

Ac 

Ac 

44 

Ac 
Ac 
s.l 
Ac 

R 

Ac 

C.1.1 

R 

Ac 



2-39JJ = 2- 



39 



2-39^ = 2-3»(N-M+l) 

2-39JJ 



2-3^ 

2-19,^ 

E^ = 2-1^+2-981? 

(p+K-1)^ 
(p^-l)o 

(q-i)o 
(q-p-K)o 

(q-p-E)e 

b 

b-a 

b-a 

1 - ^c 
N-1 



(N-1), 



b:^ 

N-1 

N-1 

b:^ 
N-1 
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VIII, 11 

12 

13 

14 

15 
16 
17 
18 



(to IX, 1) 



IX,1 
2 
3 
4 
5 
6 



8 

( 
X,l 

2 

3 

4 

5 

6 



(toX,l) 



s.l 
s.l 



S 

X 



B.l 


h 


CI 


S 


B.6 




42 


h 


C.2 


S 


C.2 




IX, 6 


Sp 


42 


h 


C.2 


s 


CI 




- 


s 


q+i-1 


s 


Cl.l 


h 


CI 


s 


C2 




38 


h- 


B.6 


h- 


42 


h- 


42 


h- 


e' 


Cc 



(to IX, 1) 



s.l 
R 

Ac 

Ac 

CI 
Ac 
Ac 
C2 

Ac 
IX, 6 

Ac 

C2 

Ac 

B'.i 

Ac 

CI 

Ac 
Ac 
Ac 
Ac 
Ac 



2-39(E.l) 



X, = a + 



(q-l)c 
1o 



(q+i-l)o 
q+i-1 

(q^i)o 
(q-^i)o 



l^(b-a) 



X- + ^ 

1 N-1 



(q^)o 

(q+i-M+D^ 
(i-M+2)^ 
(i-M+D^ 
(i-M)„ 



Note, that the box VI required no coding, hence its inmediate successor 
(VII) must follow directly upon its immediate predecessor (III). 

The ordering of the boxes is I, II, IV; III, VII, VIII, IX, X; V and 
VII, IX must also be the immediate successors of V, X, respectively, and V,2 must be 
the immediate successor of IV. This necessitates the extra orders 



V.3 

X,7 



VII, 1 
IX,1 



c 
c 



and 



IV,2 V,2 C 

As indicated in Figure 10.4, e' is 0. 
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We must now assign B. 1-6, C. 1-2, 1.1, s. 1-2 their actual values, pair the 
82 orders 1,1-25, 11,1-3, 111,1-7, IV, 1-2, V,l-3, VII, 1-9, VIII, 1-18, IX, 1-8, X,l-7 to 
41 words, and then assign I,l-X,7 their actual values. We wish to do this as a continua- 
tion of the code of 10.5. We will therefore begin with the number 53. Furthermore the 
contraits of C. 1-2, C. 1.1, s. 1-2 are irrelevant like those of 48-52 there. Hence they 
may be made to coincide with these. We therefore identify them accordingly . Siunming 
all these things up, we obtain the folloving table: 



1,1-25 
11,1-3 
IV, 1-2 
111,1-7 



53 -65 


VII, 1-9 


65 '-66' 


VIII, 1-18 


67 -67' 


IX, 1-8 


68 -71 


X,1.7 



71 '-75' 
76 -84' 
85 -88' 
89 -92 



V,l-3 
B.1-6 
C,l-2 
C.1.1 
s.1-2 



92 '-93' 
94 -99 
48 -49 
50 
51-52 



Now we obtain this coded sequence: 



49 S 
- S 
48 S 
38 h- 
42 h- 

Cc 

50 
71 C 



49 . 
42 h , 

48 , 

50 h , 

49 , 
99 h-, 
42 h-, 
85 C , 
48 S , 
a 

b 

(N-l)o 

2-39 

(P-I)o 
(q-l)o 



For the sake of completeness, we restate that part of 0-52 of 10.5 
which contains all changes and all substitutable constants of the problem. 
This is 38-52: 



48 

49 

50 
51 
52 



53 


95 


, 94 h- 


69 


42 h , 


50 Sp' 


85 


54 


51 S 


, 45 


70 


48 


50 h- 


86 


55 


94 h- 


51 - 


71 
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The 



I: 1,220 |i, II: 125 \i 
K: 300 p, X: 275 [x. 



durations may be estimated as follovs: 



III: 275 n, IV: 150 n, V: 125 u, 



VII: 350 ^t, VIII: 915 [x, 
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Total: I + II + (IV or (III + V) or III) + MI + VIII + (IX + X) x M = 
maximum = (1,220 + 125 + 400 + 350 + 915 + 575 M) ^ = 
= (575 M + 3,010) |i -/ (.6 M + 3) m. 
Hence the complete interpolation procedure (10.5 plus 10.7) requires 

(.9 M2 +.2 M + 2.4) m. 

10.8 We consider next Problem 13. h. We are again looking for that 
k = 1,...,N-M+1, for which Y^' • • • 'Yk-HW-l ^^^ ^^ close to y as possible (cf, the end 
of 10.6), i.e., for which the remoter one of yj^ and Yi^^^.i lies as close to y as 
possible. 

Since Yif-jY}^ are not equidistant, we cannot find k by an arithmetical 
criterium, as in 10,7. We must proceed by trial and error, and we will now describe a 
method which achieves this particularly efficiently. 

We are then looking for a k = E \rtiich minimizes |Jj^ = Max (Y"Ylr> Yir4M-l"y^ 
in k = 1,...,N.M+1. 

Let us first note this: Y'Yk '^ Y"Yk+l> hence y-yj^ ^ yk-HVl'^ implies 
Y-Yk ^ Mk+1» 8^d therefore M|^ > Mj^+j. Oa the other hand Yk4M-l"y ^ yk-HW'^' hence 
y^^k - ^k-^f-y inqplies M^ ^ Yk4M"y' ^^ therefore M^ ^ M^+l- Hence M^ > or g Mj^+i, 
according to whether y-yj^ > or g Y^-fWy ^'^' * ^k ^ ^k-^M* 2y < or J 0. 

In order to keep the size between -1 and 1, it is best to replace 
Vk ^ Yk^M-y - 2y by 

^i> ^k = I yk "^ I yk-^M - y • 

zj^ is monotone increasing in k. Therefore zj^ < implies zj^ "^ (^ and hence \J^ > ^-^1 
for all h i k, i.e. , |i > \x > . . . > Mj^ > M^+l- Similarly z^ - ^ implies zj^ ^ and 
hence M^ < \Ji.^\ for all h ^ k, i.e., |J|^ - Mir+l = ••• ^ 4^-M - ^-M+1* f^^om these we 
may infer 

zj^ < implies E > k, 

zk ^ implies that we can choose f g k. 

Consequently we can obtain E by "bracketing" guided by the sign of z^. 

Note that zj^ can be formed for k = 1, ...,N-M only, but not for k = N-M+1. 
llie "bracketing" must begin by testing the sign of zi for k = 1, if it is +, then E = 1, 
if it is -, then we continue. Next we test the sign of zj^ for k = N-M, if it is -, then 
E = N-M+1, if it is +, then we continue. In this case we know that 1 < E ^ N-M. Put 
k = 1, kjL " N-M. Consider more generally the case where we know; that k- < E ^ k^ . 
This implies k^-k^ § 1. If k^^-k^ = 1, then k = k-; if k^^-k^ > 1, then we continue. 
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In this case we use the function [w] which denotes the largest integer ^ w. Put 

(3) k9 = [ 1 (kt H- k:) ] , 

and test the sign of z^ for k = k?, if it is + then K ^ kP, if it is - then E > k?. 
We can therefore put kj^.^ = k"^, k^+j = k^ or k^^]^ = k?, k"!"^ ' = k^ , respectively. 
This completes the inductive step from i to i+1. Sooner or later, say for i = i , 
k^ - kj^ = 1 will occur, and the process will terminate. 

Note, that kl^i'kl^i ^ |(kt - kT -f 1), i.e., kt-kj ^ 2(k24.rk7^i^l, 
i-^-> ^i-l"^i-l = 2^^i - ^P"-^- ^^^ i "" i I however k^-k^ = 1, k^^j^-k^^.^ f 1, 

hence kt^^.k^Il ^ 2. Thus k^ -^ - k- ^j ^ 2, kt ^3 - k]; .2 ^ ^» •• • » ^V^l ^ 2^°" +1. 

00 00 

However, kj^-kj "" N-M-1, therefore, 

(4) i^ g 2iog(N.M.2) + 2. 

The virtue of this ''bracketing" method is, of course, that the number 
of steps it requires is of the order of ^log N, and not, as it would be with most 
other trial and error methods, of the order of N . (Note, that N is likely to be 
large compared to M, which alone figures in the estimates of 10.5 and 10.7.) 

After Ti has been obtained, we can utilize the routine of Problem 12 
to complete the task, just as at the corresponding point of the discussion of 
Problem 13. a in 10.7. This means that we propose to use the coded sequence 0-52 
of 10.5, and that we will adjust the coded sequence which will be formed here, 
just as in 10 .7. 

Tbe situation with the constants of Problem 12 is somewhat simpler 
than it was in 10.7. Again, p, q, and x need be given values which correspond to 
the new situation, x is clearly our present y. p,...,p+M-l and q,...,q"^-l are 
the positions of the pj, .••,P[^ and xi,...,xj^ of Problem 12, i.e., the positions of 
our present qj^ » • • • » qC+JVI- 1 ^^ ^E* ' • • '^Ic-HVi-l* ^•^•' ^^^^ ^^^ p+E-1, . . . ,p+E+M-2 and 
q+E-l, . . . ,q+Ic+M-2. Hence p = p+E-1, q = q+E-1. Tbe complications in 10.7, connected 
with q, or, more precisely, with the yi7j • • • »yi74M.l> clo ^^^t arise here, since Problem 
13. b provided for storing the entire sequence yi,.-.,>^« 

In assigning letters to the various storage areas to be used, it must 
be remembered, just as at the corresponding point in 10.7, that the coded sequence 
that we are now developing is to be used in conjunction with (i .e. , as a supplement 
to) the coded sequence of 10.5. As in 10.7, we have to classify the storage areas 
required by the latter, but this classification now differs somewhat from that one 
of 10.7: We have the storage areas D-F, which are incorporated in the final enumera- 
tion (they are 38-51 of the 0-52 of 10.5); A and B (i.e., p, . . . ,p-HVl.l and q, . . . ,q+M-l), 
which will be part of our present A and B (i.e., of p, ...,p-+N-l and q, . . . ,q-^N-l, they 
will be p+E-l,...,p+C-+M.2 and q+E-l,...,q+E+M-2); and C (i.e., r, . . . ,r-+M-2) , which 
may be at any available place. Therefore, we can, in assigning letters to the 
various storage areas in the present coding, disregard those of 10.5, with the 
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exception of C, Furthermore, there will again be no need to refer here to C (of 
10.5), since the coded sequence of 10.5 assumes the area C to be irrelevantly 
occupied. We will, at any rate, think of the letters which are meant to desig- 
nate storage areas of the coded sequence of 10.5 as being primed. We can now 
assign freely unprimed letters to the storage areas of the present coding. 

Let A and B be the storage areas corresponding to the intervals from 
p to p-W-1 and from q to q^-1. In this way their positions will be A. 1,...,N 
and B. 1,...,N, where A. i and B. i correspond to p+i-1 and q+i-1 and store q- and 
y^. As in previous problems, the position of A and B need not be shown in the 
final enumeration of the coded sequence. 

The given data of the problem are p, q, M, N, y, also r of the -coded 
sequence of 10.5. M, r are already stored there (as (M-1) , tq at 38, 43); and 
y, too (it coincides with x at 45). The definitions of p, q involve E (p^p+Ic-1, 
q = q+Ic-l), and Ic originates in the machine. The form in which p is stored accord- 
ingly involves E (as p = (p"*-E-1)q at 41); while the form in which q is stored does 
not happen to involve E (as (q-p-1) = (q-p-1) at 44). Hence 41 must be left empty 
(or rather, irrelevently occupied) when our present coded sequence begins to operate, 
and it must be appropriately substituted by its operation. 44 might be used to store 
(q-p-1) from the start, but we prefer to leave it, too, empty (i.e., irrelevantly 
occupied) and to substitute it in the process. Hence the constants requiring storage 
now (apart from M, r, y which are stored in 0-52, cf. above) are p, q, N. They will 
be stored in the storage area C. (It will be convenient to store them as (p-1) , 
(q-1) , (N-1) . ) The exit- locations of the variable remote connections will also 



be accommodated in C. The induction index i need not be stored explicitly, since 
the quantities kj, k^, k^^ contain all that is needed to keep track of the progress 

of the induction. (The i = i for which the "bracketing", i.e., the induction, 
ends, is defined by k- -k^ = 1, cf. above.) k-, k^, k- are position marks, they 

will be stored as (q+k?-l)^, (q+k^-l)^, (q+kt-l)^ (i.e., (B.k?)^, (B.kp^ ^^'KK^ 

in the storage area D. Finally, k, when formed (in the form (q+E-l) , i.e., 
(B.E) ) , will be stored in D, replacing the corresponding expression of k^. 

The index i runs from 1 to i . In drawing the flow diagram, it is 
convenient to treat the two preliminary steps (the sensing of the sign of Zj^ 
for k = 1 and k = N-M) as part of the same induction, and associate them 
ideally with two values of i preceding i = 1, i.e., with i = -1 and i = 0. 
For i = -1, k^ may be defined as the value of k for which the sign of Zj^ 
is being tested (this is its role for i > 1), i.e., as 1, N-M, respectively; 
while k^, k^ may remain undefined for i = -1, 0. 

Our task is to calculate k (in the form (q+E-1)^, i.e., (B.C)^); to 
substitute p = (p+E-1) and (q-p-1) = (q-p-D into 41 and 44; and finally to 

*^o '^ O ^ ^ O ^ 

send the control, not to e, but to the beginning of 0-52, i.e. to 0. 

We can now draw the flow diagram, as shown in Figure 10.5. The 
variable remote connections CX and (3 are necessary in order to make the differ- 
ing treatments for i = -1, for i = 0, and for i ^ 1, possible (cf. above). 
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Regarding the actual coding we make these observations: Since kj^, k^ 
are undefined for i = -1, 0, therefore the contents of D.2, 3 are irrelevant for 
i = -1, (cf. the storage distributed at the middle bottom of Figure 10.5). Ihe 
forming of (B.k^) in X, with ^^ = |W (k^+kp ], is best effected by detouring over 

the quantities 2"39(q+k£-l), 2"39(q+k£-l), 2'3s(q-fk£-l), because the operation m = [s n] 
is most easily performed with the help of 2'^^^ 2"3Qn: Indeed the right-shift R carries 
the latter directly into the former. The transitions between X and 2''39^ (i^ both 
directions) are performed as discussed at the corresponding place in 10.7. 

We will use 0-52 in 10.5 just as in 10.7, and the remarks which we made 
there cm this subject apply again. 



The static coding of the boxes I -XI follows: 
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The ordering of the boxes is I, II; III, XI; IV; V, IX, X; VI; VII; 
VIII, and II, XI, IX, IX, II must also be the inmediate successors of IV, VI, VII, 
VIII, X, respectively. Ihis necessitates the extra orders 
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As indicated in Figure 10.5, e^ is 0. 

0^1, 0^2, 0^3, 3i, (32, Ps correspond to 111,1, V,l, VII, 1, IV, 1, VI, 1, 
VIII, 1, respectively. Hence in the final enumeration 111,1, V, 1, VII, 1 must have 
the same parity and IV, 1, VI, 1, VIII, 1, must have the same parity. 

We must now assign C.1-9, D.1-3, s.l their actual values, pair the 78 
orders 1,1-7, 11,1-14, 111,1-3, IV, 1-9, V,l-ll, VI, 1-6, VII, 1-3, VIII, 1-3, IX, 1-4, 
X, 1-9, XI, 1-9 to 39 words (actually two dummy orders, necessitated by the adjusting 
of the parities of V, 1 and VIII, 1 to those of 111,1 and IV, 1, respectively, increase 
this to 40) and then assign I,1-XI,9 their actual values. We wish to do this again 
as a ccntinuation of the code of 10.5. We will therefore again begin with the 
number 53. Furthermore the contents of D .1-3, s. 1 are irrelevant like those of 
48-52 there. Hence they may be made to coincide with four of these. We therefore 
identify them with 48-51 there. Summing all these things up, we obtain the 
following table: 
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Now we obtain this coded sequence: 
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80 



The state of 0-52 of 10.5 nnist be exactly as described at the end of 



The durations may be estimated as follows: 

I: 275 ^1, II: 485 n , III: 125 u , IV: 350 n , V: 425 ix , VI: 225 |i , 
VII: 125 [i , VIII: 125 \i , K: 150 t^ , X: 330 \i , XI: 350 \i. 
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Total: I + [II or II X 2 or II X (i +1)] + [m or (IV + VI) or {(IV + V + 

o 

+ (VII or VIII) X (i .1) + K X i + X X (i -1))]+ XI 

O o 

(the two first alternatives in each bracket [ ] refer to two possibilities which 
can be described by io = 1 and io = 0) 

maximum = (275 + 485 (i +1) + (775 + 125 (i -1) + 150 i + 330 (i -1) +350)m. = 



= (1,090 i + 1,430) n. 

maximum = (1,090 ^log (N.M.2) + 3610) ^ji « (1.1 ^og (N-M.2) + 3.6) m. 

Hence the complete interpolation procedure (10.5 plus 10.8) requires 

(.9 M2 . .4 M + 1.1 ^og (N-M-2) + 3) m. 

10.9 To conclude, we takeup a variant of the group of Problems 12-13, 
which illustrates the way in which minor changes in the organization of a problem can 
be effected ex post, i.e., after the problem has been coded on a different basis. 
(In this ccauiection, cf. also the discussion of 8.9.) 

The description of the function under consideration in Problems 12 and 
13. b, i. e. , the two sequences p-i , . . . ,p|^ and xi , . . . jXj^ on one hand, and qj, . . . ,q^ and 
yn,...,)^ on the other, were stored in two separate systems of consecutive memory 
locations: p, . . . ,p-+^1-l and q, ...,q-^-l on one hand, and p, ...,p-*N-l and q, ...,q^N-l 
on the other, (lliese are the two storage areas A and B of those two problems. ) For 
Problem 13. a the situation is different, since here qi,...,qjsj alone are stored (at 
the locations p, . . . ,p■^N-l, storage area A). 

Assume now, that these data are stored together with the p^, x^, or the 
q^, y^, alternating. I.e., at the locations p, . . . ,p+2M-l, or p, . . . ,p+2N-l, with p^, 
X£ at p+2i-2, p+2i-l,or q^, y^ at p+2i-2, p+2i-l. 

Since this variant cannot arise for Problem 13. a (cf. above), and since 

its discussion for Problem 13. b comprises the same for Problem 12 (because the coding 

of the former requires combining with the coding of the latter, cf. 10.8), therefore 
we will discuss it for Problem 13. b. 
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We state accordingly: 

PROBLEM 13. c. 

Same as Problem 13 in its form 13. b, but the quantities yi»««->yN and 
«llf««jq{N{ sire stored in one system of 2N consecutive memory locations: 

p, p-^-1,.. .,p+2N-l in this order: q^, yi,.--,cjyj, Y^- 

We must analyze the coded sequences 0-52 of 10.5 and 53-101 of 10.8 
which correspond to Problems 12 and 13. b, and determine what changes are necessi- 
tated by the new formulation, i.e., by Problem 13. c. 

Let us consider the code of 10.8, i.e., of Problem 13. b, first. The 
changes here are due to two causes: First, y- is now stored at p+2i-l instead of 
q+i-1; second, since the code of 10.5, i.e., of Problem 12, will also undergo changes, 
the substitutions into it v/ill change. 

The first group is most simply handled in this way: Assume that p is 
odd. Replace the position marks (B. i) , which should be (p+2i-.l) instead of 

(q+i-1) , by their half values: (i (p+2i-l)) = (Eli + i-1) . We note, for later 

use, that this has the effect that D. 3 will contain immediately before XI 

(i (p+2E-l)) instead of (B.K) . 
2 

Returning to the general question: we must chai^ge the coding of all 
those places, where such a position mark is used to obtain the corresponding y^j^. 
The values of these position marks must be doubled before they are used in this 
way. Apart from this, however, we must only see to it that q is replaced by 

E|i . This affects C, 3 (of 10.8), i.e., 95. 

The use of position marks in the above sense occurs only in II, when 
yi and y\.^ are obtained. This takes place at 11,6 and 11,9. However, the position 
marks in question originate in the substitutions 11,2 and 11,5, and it is therefore 
at these places that the changes have to apply. 

We now formulate an adequate coding to replace II. 2-4, so as to give 
11,2 and 11,5 the desired effect. Note that at this point D.l contains 

(i (p+2k9-l)) = (Eii + k?-l)^ instead of (q+k^-D^ = ^^'^i^o* ^^''''^ ^^^"^ "'"^ 

the accumulator contains (- (D+2k--l)) • 

9 *^ 1 o 
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,1.1 




L 


2 


11,6 


Sp 


3 


38 


h 


4 


38 


h 


5 


42 


h 


6 


42 


h 



Ac 


(p.2k9-l)^ 


II. 6 


p+2k?-l 


Ac 


(p+2k?-*I-2)^ 


Ac 


(p+2k?+2M.3)^ 


Ac 


(p+2k?+2M-2)^ 


Ac 


(p+2k?+2M-l)^ 



Tlius we must replace the three orders II, 2-4 by the six orders 11,1.1-6. 
Hence the third and second remarks of 8.9 apply: We replace 11,2-3 by 11,1.1-2, then 
11,4 by 

11,1.2,1 11,1.3 C 

and let 11,1.3-6 be followed by 

11,1.7 11,5 C 

Then 11,1.3-7 can be placed at the end of the entire coded sequence. We will give a 
final enumeration that expresses these changes after we have performed all the changes 
that are necessary. 

The second group must be considered in conjunction with the code of 10.5, 
i.e., of Problem 12. In the code of 10.8 only XI is involved in the operations of 
this group. In the code of 10.5 the situation is as follows: Replace the position 
marks A.i, B. i (of 10.5), v^ich should be (p+2i-2) , (p+2i-l) , i.e., (p+2E+2i.4) , 



(p+2ie+2i-3)^, instead of (p+i-D^ 
by their half values as possible. 



second expression: 



(q+i-1) , i.e., (p+Ic+i-2)^, (q+K:+i-2)^, as nearly 
Since p is odd, we will use the half value of the 

Then we must change the 



(i (p+2E+2i-3)) =(Eii + le + i-2) 
2 2 o 

coding of all those places where such positioi marks are used to obtain the corres 

ponding p- and y- : The values of these position marks must be doubled before they 

are used for a y- , and then decreased by 1 before they are used for a pj . We must 

also take care that the storage locations, which supply the modified (A.i) , (B. i) 

position marks, are properly supplied themselves. Apart from these, however, we 

1 rr4.ou iW = (£jT + j^ + £.2) 

2 



must only see to it that XI (of 10.8) produces (- (p+2k-l)) 



instead of (p+Ic-1) = (p+E+i-2) L=]^. Since p, q, or p 
pendently, it is not necessary to produce (q-p-1) in XI 



o!i=l 
q, no longer appear inde- 



The use of position marks in the above sense occurs only in III, V, 
and VII (of 10.5), when p]^, p^, and x^, x^^.^ are obtained. Ihese take place at 
111,3, V,2, and VII, 14, 15,21. However, the position marks in question originate 
in the substitutions 111,2, V, 1, and VII, 8, 9,13, and it is therefore at these 
places that the changes have to apply. 
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Tlie position marks (A. i) , (B. i) are supplied from E.2, which in 
turn is supplied by II and IX (of 10.5). II will now supply to E.2 

(i (p+2R+l))^ = (^ +r+ i-2)^|i=2 instead of (p+1)^ = (p+E)^ = (p+I^+i-2) Ji=2- 

The function of DC is actually performed by VIII, and it will maintain in E.2 

(1 (p-^-2ir+2i-l))^= (2^ + E + i-1)^ instead of (p+i(l)^ = (p+E+i-D^. If we keep 

these facts in mind when modifying III, V, VII, then we need not change II, IX 
(i.e., VIII). 

We now give adequate codings, to insert before 111,2 and before V,l, 
and also to replace VII, 7 before VII, 8, 9, and to replace VII, 10-12 before VII, 13, 
so as to give 111,2, V, 1, and VII, 8, 9, 13 the desired effect. Note that at these 
points, as we observed above, D.4 (i.e., 41, substituted from XI of 10.8) contains 



d ,rr. 



'1 /-, 



(- (p+2ir-l)) , and E.2 contains (^ (p+2ir+l)) and (1 (p+2I?+2i.l))o, respectively. 
Hence before 111,2, V, 1, and VII, 7, the accumulator contains (-• (p+2Ic-l)) , 
(i (p+2JF+2i-l)) and (i (p+2E+2i-l)) , respectively. 



111,1.1 




L 


Ac 


(p+2E-l)^ 


2 


D.5 


h- 


Ac 


(p+2E-2)° 


V.O.l 




L 


Ac 


(p+2E+2i.l)^ 


2 


D.5 


h- 


Ac 


(p+2I?+2i-2)^ 


VII, 6.1 


D.5 


h- 


Ac 


(K (p+2E+2i-l)-l)^ 


2 ■ 




L 


Ac 


(p+2E+2i-3)^ 


VII, 9.1 


E.2 




Ac 


(^ (p+2E+2i-l))^ 


2 


D.l 


h 


Ac 


(K (p+2E-2i-l)4M-l)^ 


3 


E.l 


h- 


Ac 


('A (p+2E+2i-l)+h-l)° 


4 




L 


Ac 


(p+2K+2i+2h-3) 



Thus we must insert the orders 111,1.1-2 before 111,2; insert the orders 
V, 0.1-2 before V,l; replace the order VII, 7 by the two orders VII, 6. 1-2; and replace 
the three orders VII, 10-12 by the four orders VII, 9. 1-4. Hence the third and second 
remarks of 8.9 apply again: We replace 111,1 by 

111,0 111,1.0 C I 

and let 111,1.0 coincide with 111,1 and be followed by 111,1.1-2 and then by 

111,1.3 111,2 C I 

Next we replace V,l by 

V,0 V,0.1 c I 

and let V, 0.1 be as specified above, and be followed by V, 0.2, then by V, 1, and 
finally by 
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V,l.l 



V,2 



Finally we replace VII, 7 and VII,10.12 by VII, 6.1-2 and VII,9.1-4. It is best to 
effect this as a replacement of the entire piece VII, 7-12 by VII, 6. 1-2, VII, 8-9, 
VII, 9. 1-4 (a total of six orders to be replaced by a total of eight orders). This 
means that we replace VII, 7-11 by VII, 6.1-2, VII, 8-9, VII, 9.1, then we replace VII, 12 
by 



VII, 9.1.1 



VII, 9.2 



and let VII, 9. 2-4 be followed by 
VII, 9.2. 5 VII, 13 C I 

Then 111,1.0-3; V, 0.1-2, V,l, V,l.l; VII, 9. 2-5 can be placed at the end of the entire 
coded sequence. We will give a final enumeration that expresses these changes after 
we have performed all the changes that are necessary. 

Returning to the second group in the coding of 10.8, we note that it 
affects only XI (of 10,8). We saw above that XI must produce (- (p+2K-i)) and 
substitute it into 41, and that this is all that has to be done there. However, we 
noted before that D. 3 before XI contains precisely (i (p+2lc-l)) . Hence XI, 1-9 
(i.e., XI in its entirety) may be replaced by 



XI. 0.1 


D.3 




2 


41 


S 


3 


e' 


c 



Ac 
41 



(j (p+2IJ-l))^ 
(1 (p+2E-l))° 



Hence we replace XI, 1-3 by XI, 0.1-3, and since operations of the code of 10.8 end at 
this point, XI, 4-9 may be left empty. 



^D-1^ 



D+1 



To conclude, we observe that C. 3 must contain (El_.) = (E — - D^ 

2 2 o 

instead of (q-1) , and that neither p nor q are needed in any other form, so that 
C. 9 may be left empty. 



We are now in possession of a complete list of all changes of both 
codes 0-52 of 10.5 and 53-101 of 10.8. We tabulate the omissions and modifications 
first, and the additiais afterwards. 

Qui ss ions and Modifications: 





From 


To 


(10.8) : 


II.2-3 


11,1.1-2 




II. 4 


11,1.2.1 


(10.5) 


111,1 


111,0 




v,i 


V,0 




VII, 7-8 


VII, 6.1-2 




VII, 9-10 


VII, 8-9 




VII, 11 


VII, 9.1 




VII, 12 


VII, 9.1.1 
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Clhiissions and Modifications (cont. ): 



From 



To 



(10.8) 



XI, 1-3 

XI, 4-9 

C.3 
C.9 



XI,0.1-3 
Empty 

(Eli) 
2 o 

Empty 



Additions: 




(10.8) : 


11,1.3-7 


(10.5) : 


111,1.0-3 


(10.5) : 


V, 0.1-2 

v.i 

V.l.l 


(10.5) : 


VII, 9.2-5 



We can use five of the six empty fields XI, 4-9 (from 10.8) to accommodate 
11,1.3-7 (from 10.8) -- say XI, 4-8. C.9 is the last word of the code: 101, hence we 
can place 111,1.0-3; V,0.1-2, V,l, V.l.l, VII,9.2-5 (from 10.5) in a sequence that 
begins there. They will then occupy these positions: 



111,1.0-3 

V. 0.1.2 



101-102' 
103-103' 



V.l 
V,l.l 



104 
104' 



VII, 2-5 



105-106' 



(All these are from 10.5.) 



Now we can formulate the final form of all changes in the coded sequences 
0-52 and 53-101: 



5 
10 
14 
15 
16 
17 
57 



101 c , 




103 C , 






42 h- 


L, 


18 Sp' 


21 Sp'. 


49 


105 C , 




L , 


59 Sp 



58 


66 C 




95 


(^>o 


65 


50 


, 41 S 


101 


41 ^ , L 


66 


e'C 


, 38 h 


102 


42 h-, 5 C 


67 


38 h 


42 h 


103 


L , 42 h- 


68 


42 h 


, 58 C 


104 


10 Sp', 10 c 


69 


-- 




105 
106 


38 h , 48 h- 
L , 17 C 



The following durations are affected: 



(10.5) 



(10.8) 



III: 


+ 130 n 


V: 


+ 130 n 


VII: 


+ 130 n 


II: 


+ 180 M. 


K: 


- 225 p. 
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Total: (10.5) : (130 + 130 (M-1) + 130 MiMzii ) ^ = 

= (65 M= + 65 M) ti» (.07 M^ + .07 M) m. 
Total: (10.8) : (180 (i + 1) - 225 i ) ti = (- 45 i + 180) n 

O 

maximum SS .2 m. 
Hence there is no relevant change in the estimate at the end of 10.8. 
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11.0 CODING OF SOME COMBINATORIAL (SORTING) PROBLEMS 



11.1 In this chapter we consider problems of a combinatorial, and 
not analytical character. This means, that the properly calculational (arith- 
metical) parts of the procedure will be very simple (indeed almost absent), and 
the essential operations will be of a logical character. Such problems are of 
a not inconsiderable practical importance, since they include the category usu- 
ally referred to as sorting problems* They are, furthermore, of a conceptual 
interest, for the following reason. 

Any computing machine is organized around three vital organs: The 
memory, the logical control, and the arithmetical organ. The last mentioned 
organ is an adder- subtractor-multiplier-divider, and the multiplier is usually 
that aspect which primarily controls its speed (cf. the discussions of Part I.) 
The efficiency of the machine on all analytical problems is clearly dependent 
to a large extent on the efficiency of the arithmetical organ, and thus primarily 
of the multiplier. Now the non-analytical, combinatorial problems, to which we 
have referred, avoid the multiplier almost completely. (Adding and subtracting 
is hardly avoidable even in purely logical procedures, since the operations are 
needed to construct the position marks of memory locations and to effect size 
comparisons of numbers, by which we have to express our logical alternatives.) 
Consequently these problems provide tests for the efficiency of the non-arith- 
metical parts of the machine: The memory and the logical control. 

11.2 We will consider two typical sorting problems: Those of 
meshing and of sorting proper. Ihere are, of course, many others, and many 
variants for each problem, but these two should suffice to illustrate the main 
methodical principles. 

In order to formulate our two problems, we need the definitions which 
follow. 

We operate with sequences of complexes, A complex X = (x; u^ , . . . ,u ) 
consists of p+1 numbers: Ihe principal number x and the subsidiary numbers 
u ,...,u . p is the order of the complex. A sequence S= (X^^ , . . . ,X^° '^ 
consists of n complexes X^^^ = (x^^^; u^^^ ^ . . . ,u^ ^ M, n is the length of the 
sequence. Throughout what follows we will never vary p, the (arbitrary but) fixed 
order of all complexes to be considered. We will, however, deal with sequences 
of various lengths. 

A sequence S = (X^ ^ \ . . . ,X^° ') is a monotone if the principal elements 
x^*^ of its complexes X^ * ^ = (x^^^ u|^\...,u^^M form a monotone, non decreas- 
ing sequence: 

(1) x^l) g x<2) g ... g x<"' 

A sequence T - (Y'^ \ . . . ,Y^" is a permutation of another sequence S = (X , 
...,X^'*Mof the same length, if for a suitable permutation i— ^i' of the 
i = 1,... ,n 

(2) Y^i) = X(i'> (i = l,...,n). 
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Given two sequences S = (X^^ \ . . . ,X^°>), T - (Y^^ \ . . .,Y^»>) (of not necessarily 
equal lengths n, m), their sum is the sequence [S, T] = (X^^ \ . , . ,X^° \ Y^^ \ . . . ,Y^" ^) 
(of length n + m). 

We can now state our two problems: 

PROBLEM i4. (Meshing) 

Two monotone sequences S, T, of lengths n, m, respectively, are stored 
at two systems of n(p+l), m(p+l) consecutive memory locations respectively: 
s, s+1, . . . , s"hi(p+l)-l and t, t+1,..., t-hn(p+l)-l. An accomodation for a sequence 
of length n + m, coisisting of (n + m)(p + 1) consecutive memory locaticns, is 
available: r, r+1,..., r+(n"hii)(p+l)-l. Ihe constants of the problem are p, n, 
m, s, t, r, to be stored at six given memory locaticns. It is desired to find 
a monotone permutation R of the sum [S, T] , and to place it at the locations 
r, r+1,..., r+(n^)(p+l)-l. 

PROBLEM 15. (Sorting) 

A sequence 5 (subject to no requirements of monotony whatever) is stored 
at n(p+l) consecutive memory locaticais: s, s+1,..., s-hn(p+l)-l. Ihe constants of 
the problem are s, n, p, to be stored at three given memory locaticms. It is 
desired to find a monotaie permutation S* of S, and to replace S by it. 

11.3 We consider first Problem 14, which is the simpler one, and whose 
solution can be c cnveniently used in solving Prcblem 15. 



u 



We have S= (X<1\..., X^'^^), T= (Y^^^,..., Y<»>) withX^^^ = (x ^^ ^ 
uji^..., u^i^), Y^J^ = (y^J^ v{J\..., v^J^), and we wish to form /? = (Z ^^\..., 
2(n+m)) ^it£ z^k) = (z^k). Wj^^\... ,w^^^),^which is a monotone permutation of [S, T], 

Forming /? can be viewed as an inductive process: If Z , . . . ,Z for 
a k=l,...,n+m, have already been formed, the inductive step c aisists of forming 
Z . It is preferable to allow k = n-hn+l,too, in the sense that when this value 
of k is reached, the process is terminated. Assume, therefore, that we have forned 
Z ,...,Z *^"^^ for a k = 1, . .. ,n-hn+l. It is clear, that we may assume in addition, 
that these Z^^\...,Z^^-^^ are the X^^ \ . .. ,X^ H^ together with the Y^^ \ ... ,Y^ 4^, 
for two suitable i , j with 

(1) ^k "^ Jk "" ^"^ ^k "" ^»---' ^» Jk ^ 1»---' ^' 

Then we have the follwing possibilities: 

(a) i. = n, j. = m : By (1) k = n+m+1, hence, as we <:i)served above, 
R is fully formed, and nothing remains to be done. 

(b) Not i =n, j' =m: By (1) k ^ n -^-m, hence, as we observed above, 
R is not fully formed, and Z^^' must be formed next. Clearly, we 
may choose Z^^^ as one of X^^"^^^ and Y^Jk"^^\ we must decide which. 
We introduce a quantity 0). , so that U), =0 for Z^*"^ = X^ ^"^ ^^ and 

03, =1 {otZ^^^-Y^^^-^'K 
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03^ - implies ij^+j = i^ "^ 1| j^+i = Jk^ ^k " ^ implies i^+j = i^^, 
j\+l ~ 3^ ^1. ^^ is determined according to the following rules: 

(b.l) i^ = n, hence j^^ t m: 5 is exhausted, hence is)^ = 1. 

(t>-2) jj^ = m, hence i^^ + n: T is exhausted, hence OJ = 0. 

(b.3) i^ 4= n, jj^ 4= m: Neither S nor T is exhausted. We must 
distinguish further: 



(b.3.1) x^^k^^^ < y^Jk"^^^ Necessarily b)^ = 0. 

(b.3. 2) x^^k"^^^ > y^Jk^^^ Necessarily 0)^^ = 1. 

(b.3. 3) x^^k"*"^^ = y^Jk'*'^^ Both W^ = and oo^ = 1 are acceptable. We 

choose (0, =0. 
k 

As pointed out above, (a), (b) describes an inductive process which runs 
over k = 1,..., n+m+1. It begins with k = 1 and i ~ j\ ~ ^ (owing to (1)), and 
then progresses from k to k+1, forming the necessary i. , j, and (0 as it goes along. 

It is convenient to denote the locations s and t, where the storage of the 
sequences S and T begins, by a coirmon letter and to distinguish tKem by an index: 

(2) s^ = s, s^ = t. 

Correspondingly, let A and A be the storage areas corresponding to the two inter- 
vals of locations from s to s +n(p+l)-l and from s to s + m(p+l)-l. In this way 

their positions will be A .1, , n(p+l) and A^.l,..., m(p+l), where A .a and A .b 

correspond to s ■♦'a-l and s +b-l, and store x for a = (i-l)(p+l)+l, u^^ for 
a = (i.l)(p+l)+q+l, y^J> for b = (j-l)(p+l)+l, v^J ^ for b = (j-l)(p+l)+q + 1, 
respectively. Next, let B be the storage area intended for the sequence R, corres- 
ponding to the interval of locations from r to r+(n+m)(p+l)-l. In this way its 
positions will be B.l,..., (n4-m)(p+l), where B.c corresponds to r+c-1, and is 
intended to store z^''^ for c = (k-l)(p+l)+l, w^*'^ for c = (k-l)(p+l)+q+l, 
respectively. As in several previous problems, the positions of A , A , and B 
need not be shown in the final enumeration of the coded sequence. 

Further storage capacities are needed as follows: Tlie given data 
(the constants) of the problem, s , s , r, n, m, p, will be stored in the storage 
area C. (It will be ccaivenient to store them as (s )q, (s )q, rQ, (n(p+l))Q, 
(m(p+l))Q, (p+1)q. We will also store 1q in C. ) Next the induction indices 
k, q* will have to be stored. It is convenient to store i,, j,, too, along 
with k, while it turns out to be unnecessary to form and to store W explicitly. 
These indices i, j, k, q' are all relevant as position marks, and it will prove 
convenient to store in their place (s^ +i(p+l))Q, (s^ + j(p+l) ) q, (r + (k-l)(p+l) )o 
(s^+h(pn)+q^-l)o, (r+(k-l)(p+l)+q^.l)o (i.e. (A^Kp+D+l)^, (A^ j(p+l)4-l)^, 
(B.(k-l)(p+l)+l)o, (A^.h(p+1) +qOo, (B.(k-l)(p+l)+qMo; here W = w^^, and h = i 
or j for (i) = or 1, respectively) Note; that two position marks that correspond 
to q' are being stored. The three first quantities will be stored in the storage 
area D, and the two last ones in the storage area E. 
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We can now draw the flow diagram 
coding obtains from this without a need for 



as shown in Figure 11.1. TTie actual 
further comments. 
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tatic coding of the boxes I-X follows 
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D.l 


S 


, D.l 
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Ac 
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D.2 
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D.2 
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Ac 


6 


(to 


11,1) 


D.3 
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D.3 


C.4 






(n(p+l))o 






11, 1 






D.l 




Ac 


2 






Cl 
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Ac 


3 






C.4 


h- 


Ac 


4 


(to 


111,1) 


IV, 1 


Cc 




C5 






(in(p+l))o 






111,1 






D.2 




Ac 


2 






C.2 


h- 


Ac 


3 






C.5 


h- 


Ac 


4 


(to 


V.I) 


vi.i 


Cc 




IV. 1 






D.2 




Ac 


2 






C2 


h- 


Ac 


3 






C.5 


h- 


Ac 


4 


(to 


VII, 1) 


e 


Cc 




v.i 






D.l 




Ac 


2 






V,6 


Sp 


V,6 


3 






D.2 




Ac 


4 






V.5 


Sp 


V.5 


5 






" 







s^ +j(p+l) 

h- 
s^+Kp+l) h- 

VI. 1 Cc 



(to VII, 1) 



Ac 



Ac 



(s"), 



(r)o 
(r)„ 



(sO +i(p+l))o 

(i(p+l))o 
((i-n)(p+l))g 



(s^+j(p+l))o 
(j(p+l))o 

((j-m)(p+l))o 

(j(p+l))o 
((j-m)(p+l))o 



(s^+Kp+l))^ 
s°+i(p+l) 

(s^+j(p+l))o 
s^+j(p+l) 



y,- 



j+i 



>'j+i *i+i 



h- 



-53- 



D.l (Ahl. 

z (M'.i), 
3 cai)» 



a> 



1 



(A\l)o to D.l 
(A'.a ^ 2 
(BJX> to 3 



D.I (AUCfoi-O-i-Oo 



M 



i 




3(3.(J<-0C/^fj)HX 






(AUCM-^Oo io E.I 



# 



:m: 



1 = i+i 

>i = i 



(A'.Uf>-n)+i)o to E.I 






/o 



27./ 


(^U'Cfi-O+Oo 


Z 


(A'.f%i.i) +i\o 


£./ 


(A%(/-n)--')' 




(to 

(<{=(,.-,p) 



EH 



(<i'i, p) 



;.+ 






•yzTT 



(B.i<(rf-0-*-Oo toD.3 



L^> 



te B.CK-l)(f<rH)^i 



\7 



P. 3 (B.(K'~t)(p+ihll 



I -^ 



32 



/o 



2 



• ^ 









- /or <i><l' 



IX 






(A^h ()^-f/;-^yV/;o to EJ 



^ h ^ ^' 



Figure 11.1 



B.fK-lXjoHj-hJ 






A\ i (p^i) -hi 

A*.f(p^l)'f-<fi'l 



yd) 
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C.6 


(pH)„ 




vi.i 


D.l 




2 


E.l 


S 


3 


C.6 


h 


4 


D.l 
(to VIII, 1) 


s 


VII, 1 


D.2 




2 


E.l 


s 


3 


C.6 


h 


4 


D.2 
(to VIII, 1) 


s 


VIII, 1 


D.3 




2 


VIII, 9 


Sp 


3 


E.2 


s 


4 


C.6 


h 


5 


0.3 


s 


6 


E.l 




7 


VIII,8 


Sp 


8 


^ ^ 





s"-*(p+l) ] 



r+(k-l)(p+l) S] 



(to IX,1) 



(toX,l) 



Ac 
E.l 
Ac 
D.l 

Ac 
E.l 

A« 

D.2 

Ac 
VIII, 9 

E.2 

Ac 
D.3 
Ac 
VIII,8 

Ac 



(sUi(p+l))„ 
(s''+i(p+l))o 
(s<' + (i+l)(p+l))t 
(s°+(i+l)(p+l)). 



( 



S--^j(p+l))o 
(s^+j(p+l))o 

(s^+(j+l)(p+l))o 
(si+(j+l)(p+l))o 



(r+(k-l)(p+l))o 

r+(k-l)(p+l) 

(r+(k-l)(p+l))o 

(r+k(p+l))Q 

(r+k(p+l))o 

(s<^+h(p+l))o 

s'^+h(p+l) 



B.(k-l)(p+l)+l z 



(k) 



C.7 


\ 










\X,l 


E.2 






Ac 


(r + (k-l)(p+l)+q'-l)o 


2 


C.7 


h 




Ac 


(r+(k-l)(p+l)+q')o 


3 


IX, 10 


Sp 




IX, 10 


r+(k-l)(p+l)+q' 


4 


E.2 


S 




E.2 


(r + (k-l)(p+l)+q')o 


5 


EJ 






Ac 


(s'^^li(p+l)+q'-l)o 


6 


C.7 


h 




Ac 


(s<^4h(p+l)+q')o 


7 


IX, 9 


Sp 




IX, 9 


s'^ +h(p+l)+q' 


8 
9 


E.1 


s 




E.l 


(s'*'+h(p+l)+q')o 


s'^+h 


(p+1) +q' 


] 


Ac 


q 


10 


-. 


S 










r+(k-l)(p+l)+q' 


s] 


B.(k-1)( 


p+i)+q'+i *,<■;> 
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.1 


E.2 




2 


D.3 


h- 


3 


C.7 


h 


4 


11. 1 


Cc 



Ac 


(r+(k-l)(p+l>+q')(, 


Ac 


(q'-p-Do 


Ac 


(q'-p)o 



(to IX.l) 

The ordering of the boxes is I, II, III, V, VII, VIII, IX, X, and VII, 
VIII, IX must also be the iimiediate successors of IV, VI, X, respectively. This 
necessitates the extra orders 



IV, 5 


VII, 1 


c 


VI.5 


VIII, 1 


c 


X,5 


IX, 1 


c 



We must now assign C. 1-7, D. 1-3, E.1-2 their actual values, pair the 59 
orders 1,1-6, 11,1-4, 111,1-4, IV,l-5, V,l-7, VI, 1-5, VII, 1-4, VIII, 1-9, IX, 1-10, 
X,l-5 to 30 words and then assign I,l-X, 5 their actual values. Hiese are expressed 
in this table: 



1.1-6 

11,1-4 

111,1-4 

V,l-7 



-2' 
3 -4' 
5 -6' 
7 -10 



VII, 1-4 

VIII, 1-9 
IX, 1-10 
X,l-5 



10'-12 
12 '-16' 
17 -21' 
22 -24 



IV,1.5 
VI, 1-5 

C.1-7 
D.1-3 
E.1-2 



24 '-26' 

27 -29 
30 -36 
37 -39 
40 -41 



Now we obtain this coded sequence: 






30 , 


37 S 


1 


31 , 


38 S 


2 


32 , 


39 S 


3 


37 


, 30 h- 


4 


33 h- 


, 24 Cc' 


5 


38 


, 31 h- 


6 


34 h- 


, 27 Cc 


7 


37 


, 9 Sp' 


8 


38 


, 9Sp 


9 


- 


, - h- 


10 


27 Cc 


, 38 


11 


40 S 


, 35 h 


12 


38 S 


, 39 


13 


16 Sp' 


, 41 S 



14 


35 h 


, 39 3 


15 


40 


, 16 Sp 


16 


. 


, - s 


17 


41 


, 36 h 


18 


21 Sp' 


, 41 S 


19 


40 


, 36 h 


20 


21 Sp 


, 40 S 


21 


. 


, - s 


22 


41 


, 39 h- 


23 


36 h 


, 3 Cc 


24 


17 C 


, 38 


25 


31 h- 


, 34 h- 


26 


e Cc 


, 10 Cc' 


27 


37 


, 40 S 



28 


35 h , 3' 


29 


12 C, 


30 


is'). 


31 


(s^)o 


32 


'^o 


33 


(n(p+l))o 


34 


(m(p+l))o 


35 


(pH-l)o 


36 


lo 


37 




38 


- 


39 


- 


40 


. 


41 


- 



The durations may be estimated as follows: 

I: 225 |i, II: 150 ^, III: 150 [i, IV: 200 \i, V: 275 |i, VI: 200 n, VII: 150 li, 
VIII: 350 n, IX: 375 |i, X: 200 [i. 

„ , ^ /II +(III or (III + V) or IV) + (VI or VII) +1 . , . 
Total: I + (^ VIII , (IX + X) X p / " ^^^ ""^ " (" ^ ^^^ 
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maximum = (225 + (150 +425 +200 +350 +575 Xp) x (n + m) + 350) p. = 
= ((575 p + 1,125) (n + m) + 575) n c^ 
^.6 ((p.+ 2) (n + m) + 1) m 

Tliis result will be analyzed further in 11.5. 

11.4 We pass now to Problem 15, which will be solved, as indicated 
above, with the help of the solution of Problem 14. 

We haveS = {X^]\. ... ^X^-]) withX^^^ = (x^^^ u{i \ . . . ,u^i >), and we 
wish to form an S* = (X^^ ' \ . . . ,X^°' ^), which is a monotone permutation of 5. 

This can be achieved in the following manner: If a certain permutation 
S*® = (X ,...,X° )of5 has been found, which may not yet be fully monotone, 
but in which at least two consecutive intervals X^^' \ .. . ,X^*^"^^" ^ ' ^ and 
X +* ,o..,X +^+J" are (each one separately) monotone, then we can mesh 
these (in the sense of Problem 14), and thereby obtain a new permutation 
5*oo = X^^'°°\...,X^^'°°^) of S*^ i.e. of S, for which the whole interval 

/,»00\ I \rJ-' J.' "l ' ^^ ) ' 

X , . . . ,X "^^"^J * is monotone. Thus we can convert two monotone 
intervals of the lengths i, j into one of the length i + j. The original 5 
is made up of monotone intervals of length 1, since each X^^ (separately) may 
be viewed as such an interval. Hence we can, beginning with these monotcxie 
intervals of length 1, build up monotone intervals of successively increasing 
length, until the maximum length n is reached, i.e. the permutaticm of S at 
hand is monotone as a whole. 

This qualitative description can be formalized to an inductive process. 
Hiis process goes through the successive stages v = 0, 1, 2,... . In the stage 

V a permutation S*^ =(X^ ,...,X^° of 5 is at hand, which is made up of 
consecutive intervals of length 2^, which are (each one separately) monotone. 
(Since n may not be divisible by v, the length of the last interval has to be 
the remainder r of the division of n by 2^, r < 2^. If n is divisible by 2^, 
we should put r - q and rule this interval to be absent, but we may just as 

well put r^ = 2^, and still rule it to be the last interval.) Now a sequence of 
steps of the nature described above will produce a new permutation S*^ of <S*^, 
i.e. of 5, in which the monotone intervals of length 2^ in S*^ are paiiwise 
meshed to monotcne intervals of length 2^"^^ in S*^"*"^. (The end-effect for v + 1 
can be described in the fame terms as above for v.) Ihis inductive process 
begins vith v = and S*^ = 5. It ends when a v with 2^ § n has been reached, 
the corresponding S*^ is the desired (monotcne) S*. (Note, that at this point 
there is only one interval: It is the final interval of length r^ ^ 2^ referred 
to above, in this case with r^= n . ) Using the function <z> , which denotes the 
smallest integer ^ z, we may therefore say that this process terminates with 

V = <^log n>. 

The induction over v = 0, 1, 2,..., <^ log n> is, however, only the 
primary induction. For each v we have<(~^)^ intervals of length 2^ (cf. above), 
these form<^i<^^^= Kr^^ pairs, which have to be meshed to <^~~y^ intervals 
of length 2^"''^ (cf. above), in order to effect the inductive step from v to v + 1. 
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If we enumerate these interval pairs by an index w = 1,..., <(-f^^ , then it becomes 
clear that we are dealing with a secondary induction, of which w is the index. 

We can now state the general inductive step with complete rigor: 

Consider a v = 0, 1, 2,..., and assume that the permutation 

*^ ~ vX , . . . ,X ° ) has already been formed. We effect a secondary induction 
over a w = 0, 1, 2,... . Assume that a certain w has already been reached. Then 
2"^ w ^ n. We now have to mesh two intervals of the length n^, m^, i.e. X^^'°^ with 
i = 2^+lw + l,..., 2^+iw+n^ and X^^'^^ with j + 2-+lw + n^ + 1, . . . , 2^+lw+n^+m^. 
Ordinarily %''m^ = 2^, i.e. this is the case when both these intervals can be 
accomodated in the interval l,...,n. Thus the condition fcr this case is 
2^+ (w+1) i n. Otherwise, i.e. for 2^+Uw+l) > n, we have two possibilities: 
First: n^ - 2^ if this interval can be accomodated in the interval l,...,n, i.e. 
if 2^+ w+2^ ^ n. m^ is then chosen so as to exhaust what is left of the interval 
l,...,n, i.e. m^ = n - (2^ w + 2^). Hius the condition for this case is 

2^+1 (w+1) > n ^ 2^+^w + 2\ We may include 2^+^(w+l) = n, too, since the formulae 
of the present case coincide then with those of the previous one. Second: Other- 
wise, i.e. for 2^"^ w + 2^ > n, the second interval is missing m^ = 0. % is then 
chosen so as to exhaust what is left of the interval l,...,n, i.e. n^ = n - 2^"^^w. 
Thus the condition for this case is 2^-^^w + 2^ > n (and, of course, g 2^+^w). We 

may include 2^"*" w + 2^ = n, too, since the formulae of the present case coincide 
then with those of the previous one. These rules can be sunmed up as follows: 

[%= 2v , m^=2^ 

for 2^+Mw+l) ^ n, 

(3) < % ^ ^' ' % = ^ - (2^^'w+2M 

^ for 2^-*-^(w+l) § n ^ 2^-^lw+2^, 

n^ = n-2^+^w s %= 
^ for 2^+^w + 2^ ^ n (^ 2^+lw). 

j£ 2^+ (w+1) < n, then the (secondary) induction over w continues to w + 1; if 
2^+ (w+1) ^ n, then the induction over w stqps. If the induction over w stops 
with a w > 0, then the (primary) induction over v continues to v + 1; if the 
induction over w stops with w = 0, then the induction over v stops. At any rate 
these meshings form the permutation S*^"^ . If the induction over v stops with 
a certain v, then its 5*^"^ is the desired 5*. 

There remains, finally, the necessity to specify the locations of S*^ 
and S*^+ , which control the meshing processes that lead from 5*^ to S*^'^ . 
Assume that 5*^ occupies the interval of locations from a ^ to a ^ +n(p+l)-l, 
and 5*^+^ similarly the interval of locations from a^^"^^^ to a^^"^^^ + n(p+l)-l. 
Then the meshing process which corresponds to a given w meshes the sequences that 
occupy the intervals from a^^^ +2^'^^w(p+l) to'a^^^ + (2^*^V +n^)(p+l)-l and from 
^(v) +(2v4.1^+jj^)(p+i) to a^^^ +(2^^^w+n^+m^)(p+l)-l, and places the resulting 
sequence into the interval from a^^+^^ +2^+^w(p+l) to a^^"^^^ + (2^+^w + n^+m^)(p+l)-l. 
The meshing process is that one of Problem 14. 
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We distinguish the constants of that problem by bars. The above specifications 
of the locations that are involved can then be expressed as follows: 

fs^ =a^^U 2^+lw (p+1), 

i^ =a^-> +(2^+lw+n^)(p+l), 



(4) 



T = a^^+^^ + 2^+lw (p+1), 

p = p 

Clearly the interval from a^""^ to a^''^ +n(p+l)-l and the interval from a^"""*"^ ^ to 
a ^^ +n(p+l)-l must have no elements in common. 

We need such an a^''^ i.e. an interval from a^'^^ to a^""^ +n(p+l)-l, 
for every v for which S*"" is being formed. It is, however, sufficient to have 
two such intervals which have no elements in common, say from a to a+n(p+l)-l 
and from b to b +n(p+l)-l. We can then let a^^^ alternate between the two values 
a and b, i.e. we define 



(5) a^'^ 



J- a for V even 
= b for V odd 



Hence a^^^ = a, but S* = S, hence the statement of Problem 15 requires a^^^ 
Consequently 



(6) 



s - a. 



Thus a is the s of Problem 14, and b is any such location such that the interval 
from a to a +n(p+l)-l and the interval from b to b+n(p+l)-l have no elements in 
cormion. 

We want the, final »S* to occupy the location of the original S, i.e. 
we want a^^+^^ =a^^^ =a for the final v with S*^+^ =S*. I.e. this v+1 should 
be even. We saw further above, that the induction over v can stop when the 
induction over w stops with w=0, i.e. when 2^+ ^n. This v+1 may be odd; let 
us therefore agree to perform one more, seemingly superfluous, step from v to 
v+1. I.e. we will terminate the induction over v when 

(7) 2^+1 I n and v+1 even, i.e. a^^+^^ = a. 

Hence this v + 1 is the smallest even X with 2 g n, i.e. 2X' for the smallest 
integer V with 2^^^ 2 n; i.e. 4^^ 2 n. This means that 

(8) v+1 = 2 <('*log n^ . 
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In the coded sequence that we will develop for our problem, the coded 
sequence that solves Problem 14 will occur as a part, as we have already pointed 
out. This is the coded sequence 0-43 in 11.3. We propose to use it now, and to 
adjust the coded sequence that we are forming accordingly. 

This is analogous to what we did in 10.6 and 10.7 for the Problems 13. a 
and 13. b: TTiere the coded sequence formed in 10.5 for Problem 12 was used as part 
of the coded sequences of the two* former problems. There is, however, this diffe- 
rence: In 10.6, 10.7 the subsidiary sequence (of Problem 12) was attached to the 
end of the main sequences (of Problems 13. a, b), while now the subsidiary sequence 
(of Problenj 14) will occur in the interior of the main sequence (of Problem 15) — 
indeed it is a part of the inductive step in the double induction over v, w. The 
constants of Problem 14 have already been dealt with: They mtist be substituted 
according to (4) above. 

In assigning letters to the various storage areas to be used, it must 
be remembered, just as at the corresponding points in 10.6, 10.7, that the coded 
sequence that we are now developing is to be used in ccnjunction with (i.e. as 
an extension of) the coded sequence of 11.3. It is therefore again necessary 
to classify the storage areas required by the latter: We have the storage areas 
C-E, which are incorporated in the final enumeration (they are 30-41 of the 0-41 
of 11.1); and A^J^l^andB (i.e. i°, . . . , i^ + n (p+l)-l; s \ . . . , i^ + iS (p+l)-l 
and r,..., r + (n-h5)(p+l)-l), which will be part of our present A, B (i.e. of 
a,..., a+n(p+l)-l and b, . . . , b +n(p+l)-l, they will be a^^^ +2^+Mp+l), . . . ,a^^^ + 
+ (2^^^w +n^)(p+l).l; a<^^ + (2^+^ +n^)(p+l), . . . , a^^ ^ +(2^+lw +n^+ m^)(p+l)-l, 
and a^^+1^ +2-+Mp+l),..., a^^+^ ^+(2^+lw +n^ +m^)(p+l).l, cf. (4) above). 
Tlierefore we can, in assigning letters to the various storage areas in the present 
coding, disregard those of 11.3. 

Let A and B be the storage areas corresponding to the intervals from a 
to a+n(p+l)-l and from b to b+n(p+l)-l. In this way their positions will be 
A.l,...,n(p+1) and B. 1, . . . ,n(p+l) , where A.i and B.i correspond to a+i-1 and b+i-1. 
A will store S at the beginning and S* at the end of the procedure. B may be at 
any available place, and we assume it to be irrelevantly occupied. In the course 
of the procedure A and B are the alternate values of A^^ (for a^^^ = a and b, 
i.e. for V even and odd, respectively), and S is stored at A while S ^ or 
S^^"*" is being formed. All these storages are arranged like the storage of 5 
at A in 11.3. 

The given data of the problem, a, b, n, p, will be stored in* the 
storage area C, with the exception of p. (It will be convenient to store them 
as aQ, bQ, n(p+l)Q. A will also be stored in C. ) p is also needed in the 
coded sequence of 11.3, and it is the only one of the constants of that problem 
which is independent of the induction indices v, w. (Cf . (4) above. ) Hence its 
storage in that sequence, at 35 (as (p'^D^) is adequate for our present purpose. 
The other constants of that problem depend on v, w (cf. above), hence their 
locations, 30-34, must be left empty (or, rather, irrelevantly occupied), when 
our present coded sequence begins to operate, and they must be appropriately 
substituted by its operation. 
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The induction index v will be stored in D. (It will be convenient to 
form (2^"''^(p+1))q, and to store (a^^^ )q, (a^^*^^^ )q along with it.) TTie induction 
index w is part of the expressions which have to be placed at 30 and 32 ((s )q 
and (?)o, i.e. by (4) (a^^^ +2^+^w (p+D) and (a^^"^^^ +2^+*w (p+l))o), hence it 
requires no other storage. 



We can now draw the flow diagram, as shown in Figure 11.2. The actual 
coding obtains from this quite directly, in one instance (at the beginning of 
box VI) the accumulator is used as storage between two boxes. 

The static coding of the boxes I -XIV follows. 



Cl 






«0 




2 






bo 




1,1 






35 




2 






D.l 


S 


3 






Cl 




4 






D.2 


S 


5 






C.2 




6 






D.3 


S 


( 


to 


11,1) 






11,1 






D.2 




2 






30 


S 


3 






D."3 




4 






32 


s 




(to 


111,1) 




111,1 






30 




2 






D.2 


h- 


3 






D.l 


h 


C.3 






(n(p+l))o 




111,4 






C.3 


h- 


5 


(to 


v.i) 


IV, 1 


Cc 


IV, 1 






C.3 




2 






D.2 


h 


3 






31 


s 


4 






30 


h- 


5 






33 


s 


C.4 











IV, 6 






C.4 




7 






34 


s 



Ac 


(p+l)o 


D.l 


(p+l)„ 


Ac 


^0 


D.2 


«o 


Ac 


bo 


D.3 


bo 


Ac 


(a<'>) 


30 


(a<-)) 


Ac 


(a<'-^l 


32 


(a<'-^ 



(to IX, 1) 



Ac 
Ac 
Ac 

Ac 



Ac 
Ac 
31 
Ac 
33 

Ac 
34 



(a<'>+2^+^w(p+l))o 
(2'+V (p+1)),, 

((2^+lw + 2M(p+l))o 
((2'+lw+2'-n)(p+l))o 



(n(p+l))o 
(a<^'+n(p+l))o 
(a<^>+n(p+l))o 
((n-2'+^)(p+l))o 

((n-2^+lw)(p+l))o 
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c 


'1 

I 




3 CS./?o 




z cA^i/X 

3 (A'^^'Ul 




3Z (A'^'iOo 






















IT 












Cr*-')o to D.1 

(B.I)o to 3 


,- 


# 


3 


i 3 




4 


1 


_i /-,_- 1 


d- 


JP 1 w /■ 1 












iik: 


,3 


^ 






# 












i:? 


-\ o~^w 1 








r 






' 7 


# 




-I 






__J iy^i..jw \j 


II 


>■ ^ ► . , 


_^v_, . 








3 


"1 

1 vv •+/->• 


^ "^ 








i 






« 


.8 






1 


5- 


M 






(e) 






Z'^'w-^-Z^'-r^ 




+ i 


. B i 


l~ 


\ > 








•TT' 


— ^ 


^ 










^— ^ 




jt 1 




CA"'in(pi-i}H)a-to 3» 
O io3^ 




m 




A.h/^ri\ 


ft 








I e 




3 


1 










«l 










?v+i 


->^ 






1 


31 




7 






+> 








2''"/W<-')-' 


? 


3: 


c ^ J^- 






- > 


f c: yf 
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V.I 
2 
3 
4 
5 
(to 


D.l 
33 
30 
D.l 
31 
VI, 1) 


S 

h 
S 


Ac 
33 
Ac 
Ac 
31 


(2Mp+l))o 
(2Mp+l))„ 

(a<''>+2'+lw (p+l))o 
(a<'''+G^+'w + 2^)(p+l))o 
(a<->+(2T+l^+2')(p+l)) 


vi.i 


D.2 


h- 


Ac 


((2-+1w+2-)(p+1))q 


2 


D.l 


h 


Ac 


(2-+i(w+l)(p+l))„ 


3 


C.3 


h- 


Ac 


((2^+l(w+l).n)(p+l))o 


4 


VII, 1 


Cc 






(to 


VIII, 1) 








VII, 1 


C.3 




Ac 


(n(p+l)), 


2 


30 


h- 


Ac 


((n-2'+M(p+l)-a<'>)o 


3 


D.2 


h 


Ac 


((n-2-+M(p+l))o 


4 


D.l 


h- 


Ac 


((n-(2'+lw+2^))(p+l))o 


5 


34 


S 


34 


((n-(2^+lw+2-))(p+l))o 


(to 


IX,1) 








VIII, 1 


D.l 




Ac 


(2-(p+l))j, 


2 


34 


s 


34 


(2np+l))o 


(to 


IX, 1) 








IX 


(Problem 14: 


0-29.) 






(to 


X,l) 








X,l 


30 




Ac 


(a<')+2-+lw(p+l))Q 


2 


D.2 


h- 


Ac 


(2'+lw(p+l))o 


3 


D.l 


h 


Ac 


((2»+lw + 2-)(p+l))g 


4 


D.l 


h 


Ac 


(2'+l(w+l)(p+l))o 


5 


C.3 


h- 


Ac 


((2'+nw+l)-n)(p+l))o 


6 


XI, 1 


Cc 






(to 


XIII, 1) 








XI.l 


D.l 




Ac 


(2np+l))„ 


2 


D.l 


h 


Ac 


(2-+l(p+l))o 


3 


C.3 


h- 


Ac 


((2-+l.n)(p+l))o 


4 


XII, 1 


Cc 






(to 


XIV, 1) 








XII, 1 


CI 




Ac 


«0 


2 


D.3 


h- 


Ac 


(a-a'^+i>)o 


3 


s.l 


S 


s.l 


(a-a^^+^))„ 


4 


s.l 


-M 


Ac 


-!(a.a('^«>)4 


5 


e 


Cc 






(to XIV, 1) 
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XIII, 1 


30 


2 


D.l 


3 


D.l 


4 


30 


5 


32 


6 


D.l 


7 


D.l 


8 


32 


(to 
XIV, 1 


III.l) 
D.l 


2 


D.l 


3 


D.l 


4 


D.2 


5 


s.l 


6 


D.3 


7 


D.2 


8 


s.l 


9 


D.3 



h 
h 
s 

h 
h 
s 



h 

s 

s 
s 
s 



(to 11,1) 



Ac 

Ac 
Ac 
30 
Ac 
Ac 
Ac 
32 

Ac 

Ac 

Ac 

Ac 

s.l 

Ac 

D.2 

Ac 

D.3 



(a<^>+2^+lw (p+D), 






(a*'' +(2-+V+2')(p+l))g 



a-- +2^+l(w+l)(p+l))o 

„(t) 



(a^'"+2^+Mw+l)(p+l))o 
(a(v+l)+2v+iw (p+l))o 
(3(v+i)+(2-+iw+2'')(p+l))„ 
(3(v+l)+2v+l(w+l)(p+l))o 

(3(v+l)+2T+l(^+l)(p+l))o 

(2Mp+l))o 

(2»+l(p+l))(, 

(2-+i(p+l))o 





(a(^>) 



(a<-+l)), 



(a 



(t+1) 



(a^-0 
(a<'>) 



TTie ordering of the boxes is I, II, III, V, VI, VIII, K, X, XIII; 
XI, XIV; IV; VII; XII, and IX, IX, XIV, III, II must also be the iirmediate 
successors of IV, VII, XII, XIII, XIV, respectively. In addition, IX cannot 
be placed immediately after VIII, since IX, 1 is -- but IX must nevertheless 
be the immediate successor of VIII. All this necessitates the extra orders 



IV, 8 


IX,1 


c 


VII, 6 


IX, 1 


c 


VIII, 3 


IX, 1 


c 


XII, 6 


XIV, 1 


c 


XIII, 9 


111,1 


c 


XIV, 10 


11,1 


c 



Finally in order that X be the immediate successor of IX, the e of 
0-41 (in 26) must be equal to X, 1. 

We must now assign C.1-4, D.1-3, s.l their actual values, pair the 76 
orders 1,1-6, 11,1-4, 111,1-5, IV,l-8, V,l-5, VI, 1-4, VII, 1-6, VIII, 1-3, X,1.6, 
XI, 1-4, XII, 1-6, XIII, 1-9, XIV, 1-10 to 38 words, and then assign 1,1 - XIV, 10 
their actual values. (IX is omitted, since it is contained in 0-41.) We wish 
to do this as a continuation of the code of 11.3. We will therefore begin with 
the number 42. Furthermore the contents of D. 1-3, s.l are irrelevant, like 
those of 37-41 there. We saw, in addition, that 30-35 (with the exceptions of 
30,32) are also irrelevant. We must, however, make this reservation: 30-41 
are needed while box IX operates, and during this period D. 1-3 are relevantly 
occupied, but not s.l. s.l is relevantly occupied only while boxes XII, XIV 
operate, and during this period 30-35 and 37-41 are irrelevant. Hence only 
s.l can be made to coincide with one of these, say with 41. Suirming all these 
things up, we obtain the following table: 
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1,1-6 


42 -44' 


VIII, 1-3 


54 -55 


IV, 1-8 


70 -73' 


11,1-4 


45 -46' 


X,l-6 


55 '-58 


VII. i-6 


74 -76' 


111,1-5 


47 -49 


XIII. 1-9 


58'-62' 


XII, 1-6 


77 -79' 


V,l-5 


49'-51' 


XI, 1-4 


63 -64' 


C.1-4 


80 -83 


VI, 1-4 


52 -53' 


XIV, 1-10 


55 -69' 


D.1-3 
s.l 


84 -86 
41 



Now we obtain this coded sequence: 



42 


35 


, 84 S 


57 


84 h 


82 h- 


72 


33 S 


83 


43 


80 


, 85 S 


58 


63 Cc 


30 


73 


34 S 


C 


44 


81 


86 S 


59 


84 h , 


84 h 


74 


82 


, 30 h- 


45 


85 , 


30 S 


60 


30 S , 


32 


75 


85 h 


, 84 h- 


46 


86 


32 S 


61 


84 h 


84 h 


76 


34 S 


, C 


47 


30 


85 h- 


62 


32 S 


, 47 C 


77 


80 


, 86 h- 


48 


84 h , 


82 h- 


63 


84 


84 h 


78 


41 S 


, 41 -M 


49 


70 Cc 


, 84 


64 


82 h- 


77 Cc 


79 


e Cc 


65 C 


50 


33 S 


30 


65 


84 


84 h 


80 


^0 


51 


84 h . 


31 S 


66 


84 S 


, 85 


81 


^0 


52 


85 h-, 


84 h 


67 


41 S 


86 


82 


(n(p+l))o 


53 


82 h- 


, 74 Cc 


68 


85 S 


41 


83 





54 


84 , 


34 S 


69 


86 S 


45 C 


84 


_ 


55 


OC 


, 30 


70 


82 


85 h 


85 


- 


56 


85 h- 


, 84 h 


71 


31 S 


, 30 h- 


86 


- 





26 


55 Cc', 


10 Cc' 


33 


30 


. 




34 


31 


- 




35 


32 


- 




36 



For the sake of completeness, we restate that part of 0-41 of 11.3, which 
contains all changes and all substitutable consteints of the problem. Hiis is 26 
and 30-41: 



37 
38 
39 
40 
41 



The durations may be estimated as follows: 

I: 225 u, II: 150 \x, III: 200 n, IV: 300 \i, V: 200 n, VI: 150 p., VII: 225 \i, 
VIII: 125 \i, X: 225 n, XI: 150 n, XII: 225 n, XIII: 350 n, XIV: 375 n- 
IX: The precise estimate at the end of 11.3 is 

((575 p + 1,125) (n^ + m^) + 575) \i. 

Total: Put ? = 2 <*log r^ , ii = 2 <(*log n> -<(2log i^ = or 1, 
w = ^ " iN . TKen the total is 

\2r+l/ 
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i + ^i'Jii+'Y 

T=0 ,=0 

This is majorized by 



III + (iV or (V+VI+(VII*> or VIII))) 
+ IX +X + (XI*' or XIII) 



+ XIV**' }+xii (n+1) 



(225 + '2} 150 + '2^ 

T=0 1 ,=0 



200 + 200 + 150 + 125 + (575 p + 1, 125) X 
X (n^+ni^)+ 575 + 225 + 350 

- 375 + 225 (n+D) |i . 
Since 2 (n^ + m^) = n, fi ^ 1, this is majorized by 
(225 + '1^ {(575 p + 1,125) n + l,825<-iL_N, + 425} + 75) li . 

Since 2 <Cr75fT^ ^ 2 ( -^ + 1) | n + v, this is majorized by 

((575 p + 1,125) n V + 1,825 n + 2,250 v + 300) ^ . 
We can write this in this form: 



+ 100 - 200 + 375 



where 



((575 p + 1,200 (1 + 6)) n V ) [JL^ (.6 (p + 2 (1 + 6)) n ?) m , 

V = 2<nognN , 

6 = . .06 ^ L|2 + L13 + ^ . 
V n n V 

6 is a small and slowly changing quantity: For n = 100; 1,000; 10,000 
(the last value is, of course, incompatible with the memory capacities which 
appear to be practical in the near future) we have v = 8; 10; 14 and hence 
6 = .15; .09; .05, respectively. 



* ) "*- 1 

At most once among the w in 2 . We disregard this effect for IV, which 

represents the shorter alternative. We replace VII, XI by their alternatives, 

VIII, XIII, respectively, inside the 2 . This necessitates the corrections 

VII-VIII = 100 M., XI-XIII = -200 M,, respectively, outside the ""t . (Cf. the 

w = 

third and seccHid terms from the right in the brackets { ... } of the next formula.) 



') 



r-l 



Missing once among the v in 2 • We may therefore subtract XIV - 375 |i as 

a correction outside the 2 . 

v = 



.66^ 



11.5 The meshing and sorting speeds of 11.3 and 11.4 are best expressed 
in terms of m per complex or of complexes per minute. ITie number of complexes in 
these two problems is n + m and n, respectively, hence the number of m per complex 
is c:^,6 (p+2) and ^.6 (p+2(l+6)) • 2 <^log n >;tl.2 (p+2.3) <^*log n , 
respectively. ITie number of complexes per minute obtains by dividing these numbers 
into 60,000, i.e. it is 

(9) 100.000 and 50,000 respectivelv 

To get scMne frame of reference in v/hich to evaluate these figures, one might compare 
them with the corresponding speeds on standard, electro-mechanical punch-card 
equipment. 

Meshing can be effected with the I.B.M. collator, which has a speed of 
225 cards per minute. Sorting can be effected by repeated runs through the I.B.M. 
(decimal) sorter. Tlie sorting method v/hich is then used is not based on iterated 
ineshing, and hence differs essentially from our method in 11.4. Strictly construed, 
it requires as many runs through the sorter, as there are decimal digits in the 
principal number. This number of digits is usually between 3 and 8. Since we are 
dealing in our case v/ith 40 binary digits, corresponding to 12 decimal digits, the 
use of the value 6 in such a comparison does. not seem unfair. Hie sorter has 
a normal speed of 400 cards per minute, and it can unquestionably be accelerated 
beyond this level. 50% seems to be a reasonable estimate of this potential acce- 
leration. Making the comparison on this basis, we have a speed of 7-^— = 100 

cards per minute. Hence we have these speeds for meshing and sorting in cards 
per minute: 

(10) 225 and 100, respectively. 

A card corresponds to one of our complexes, since it moves as a unit in meshing and 
in sorting. Hence the speeds of (9) and (10) are directly comparable, and they give 
these ratios: 

(11) J^ ^ -i^ and — ^, r, respectively. 

(10) ^ p+2 (p+2.3)<(nogn>' *^ 

A standard I.B.M. punch card has room for 80 decimal digits. Hiis is 
equivalent to about 270 binary digits, i.e. somewhat less than 7 of our 40 binary 
digit numbers. It is, of course, in most cases not used to full capacity. Hence 
it is best compared to a complex with p + 1 ^ 7, i.e. p =1, ...,6. For n values 
from 100 to 1,000 seem realistic (in view of the probably available "inner" 
memory capacities, and assuming that no "outer" memory is used, cf . the second 
remark at the end of this sec^tion), hence <^ log n> = 4, 5. Consequently the 
ratios of (11) become 

(11') -|^ ::^ 150 to 55 and 30 to 15, respectively. 
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In considering these figures the following facts should be kept in mind: 

First: In using an electronic machine of the type under consideration 
for sorting problems, one is using it in a way which is least favorable for its 
specific characteristics, and most favorable for those »of a mechanical punch card 
machine. Indeed, the coherence of a complex X = (x; u^,..., u ), i.e. the close 
connection between the movements of its principal number x and the subsidiary 
numbers u^,..., u , is guaranteed in a punch card machine by the physical identity 
and coherence of the punch card which it occupies, while in the electronic machine 
in question x and each u^,..., u must be transferred separately; there is no in- 
trinsic coherence between these items, and their ultimate, effective coherence 
results by synthesis fran a multitude of appropriately coordinated individual 
transfers. This situation is very different from the one which exists in properly 
mathematical problems, where multiplication is a decisive operation, with the 
result that the extreme speed of the basic electronic organs can become directly 
effective, since the electronic multiplier is just as efficiently organized and 
highly paralleled as its mechanical and electromechanical counterparts. Thus in 
properly mathematical problems the ratio of speed is of the order of the ratio of 
the, say, .1 m multiplication time of an electronic multiplier and of the, say, 
1 to 10 second multiplication time of an electromechanical multiplier, i.e., say, 
10,000 to 100,000 — while (11') above gave speed ratios of only 15 to 150. 

Second: The "inner" memory capacities of the electronic machines that 
we envisage will hardly allow of values of n (and m) in excess of about 1,000. 
When this limit is exceeded, i.e. when really large scale sorting problems arise, 
then the "outer" memory (magnetic wire or tape, or the like) must also be used. 
The "inner " memory will then handle the problem in segments of several 100, or 
possibly up to 1,000, complexes each, and these are combined by iterated passages 
to and from the "outer" memory. This requires some additional coded instructions, 
to control the transfers to and from the "outer" memory, and slows the entire 
procedure som.ewhat, since the "outer" memory is very considerably less flexible 
and less fast available than the "inner" one. Nevertheless, this slowdown is 
not very bad: We saw in 11.3 that meshing requires .6 m per number. Magnetic 
wires or tapes can certainly be run at speeds of 20,000-40,000 (binary) digits 
per second, i.e. 500-1,000 (40 binary digit words or) numbers per second. This 
means 1-2 m per number. Thus the times required for each one of these two 
phases of the matter are of the same order of magnitude. In addition to this 
the "outer" memory is likely to be of a multiple, parallel channel type. 

We will discuss this large scale sorting problem later, when we come 
to the use of the "outer" memory. It is clear, however, that it will render 
the comparison of speeds, (11'), somewhat less favorable. 

Third: We have so far emphasized the unfavorable aspects of sorting 
with an electronic machine of the type under consideration — i.e. one which is 
primarily an all-purpose, mathematical machine, not conceived particularly for 
sorting problems . Let us now point out the favorable aspects of the matter. 
These seem to be the main ones: 
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(a) All the disadvantages that we mentioned are relative ones (i.e. in 
comparison to the properly mathematical problems) — electronic sorting should 
nevertheless be faster than mechanical, punch card sorting: In our above examples 
by factors of the order of 10 to 100. 

(b) The results of electronic sorting appear in a form which is not 
exposed to the risks of the unavoidable human manipulations of punch cards: 
They are in the "inner" memory of the machine, which requires no human manipu- 
lation whatever; or in the "outer" memory, which is a connected, physically 
stable medium like magnetic wire or tape. 

(c) Hie sorting operations can be combined with and alternated with 
properly mathematical operations. This can be done entirely by coded instructions 
under the inner control of the machine's own control organs, with no need for 
human intervention and no interruption of the fully automatic operation of the 
machine. This circumstance is likely to be of great importance in statistical 
problems. It represents a fundamental departure from the characteristics of 
existing sorting devices, which are very limited in their properly mathematical 
capabilities. 



